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ABSTRACT 


The GLAS cllMte itodel is e general circulation model based on the primitive 
equations in sigma coordinates on a global domain in the presence of orography. 
The model incorporates parameteriaations of the effects of radiation* convection* 
large scale latent heat release* turbulent and boundary layer fluxes and ground 
hydrology. 

Winter and summer simulations have been carried out with this model* ana 
the resulting data compared to observations. These comparisons indicate that: 

(i) The simulated mean dynamical fields of the upper troposphere in winter and 
the lower troposphere in summer are more realistic than those of the lower 
troposphere in winter and the upper troposphere in summer. A serious problem 
is the unrealistically cold polar upper troposphere in both simulations, (ii) 

The latitude-height structure of the components of the simulated fields with 
various spatial (planetary wave* synoptic scale wave) scales and temporal (time 
mean* low frequency* medium frequency) scales shows overall agreement with the 
observations. However, the variance of low frequency planetary waves Is somewhat 
too small, (ill) Maps of local variability Indicate that the model's band-pass 
variances are well simulated except for the upper level heat and momentum 
fluxes. The local* low-pass variances are not realistic. (Iv) The simulation 
of physical processes Is generally realistic. The distribution of simulated 
evaporation* precipitation and surface stress show good agreement with observa- 
tions. The surface sensible heat flux Is somewhat too high* and the model 
shows a serious excess of supersaturation clouds In the lowest layer, (v) The 
Hadley cells In both seasons are somewhat too weak, but their location and 
depth are realistic. The model Walker circulation in the winter subtropics 
Indicates rising over an oceanic precipitation maximum and sinking over land. 
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This tschnlcsl laemorandua has been prepared as a docunentacion of Che 
design and perforaance of the GIAS Cllaate Model* Although this global general 
circulation nodel Is continually evolving, •‘he results presented here should be 
useful both as an Indication of aodeling sk<ll and as a basis for comparison 
with future model results. 
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1. Introduction 


General circulation models are mathematical tools for conducting sensi- 
tivity, predictability and observing system simulation studies. The simple 
quasi -geostrophlc, channel model with prescribed heating and hydrologic cycle 
developed by Phillips (1956) led the way for development of more complex models 
(Smagorlnsky, 1963; Mints, 1964; Manabe et al. 1965, Leith, 1564, Arakawa et 
al, 1969), which were based on the primitive equations of motion on a global 
domain, and parameterized the physical processes of radiation, large scale and 
convective condensation, boundary layer and turbulent fluxes, and ground hydro- 
logy. The primary objective of these efforts was to simulate the space-time 
structure of the stationary and transient components of the circulation. 
Although further model developments are being carried out to determine and 
correct the deficiencies of the general circulation models (GCMs), their degree 
of success in simulating the observed circulation justifies their use for con- 
ducting appropriate sensitivity, predictability and observing system simulation 
studies. 

In this paper we present the analysis of winter and summer circulations 
simulated by the GLAS (Goddard Laboratory for Atmospheric Sciences) general 
circulation climate model. The present version of the model has undergone sub- 
stantial modifications since the publication of Its earlier version seven years 
ago by Somerville et al. (1974). Results of simulations by an intermediate 
version of the model were presented by Halem et al. (1979) at the GARP/JOC 
(Joint Organizing Committee) conference on climate models. 

The kind of controlled numerical experiments that can be carried out with 
a GCM depends upon the model's ability to simulate specific features of the 
circulation. As one of our main objectives Is to conduct sensitivity studies 
to determine the role of the slowly varying boundary conditions of sea surface 
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temperature, soil moisture, sea ice and snow In determining the predictability 
of monthly and seasonal anomalies and their geographical distribution, we present 
a detailed spatial and temporal analysis of the physical state parameters and 
processes of the model atmosphere. 

The model shows a serious deficiency in simulating the mean circulation 
near the poles. The upper tropospheric temperatures near the winter poles are 
too low compared to the observations and the corresponding upper level zonal 
winds are too strong. Another deficiency is that the tropical troposphere 
warms up by 2-*3*'C during the course of the simulation. We have not yet deter- 
mined the precise cause of this deficiency. But the simulated tropospheric 
circulation in the middle latitudes and the tropics is realistic enough to 
warrant detailed analysis of its space-time variability. 

Some of the primary objectives of NASA's climate modeling program are: 

I) to conduct simulation studies which will help optimize space-observing 
systems; and 2) to determine how the global boundary conditions of sea surface 
temperature, soil moisture, sea ice and snow affect the variability and predic- 
tability of monthly and seasonal atmospheric anomalies. These studies are being 
carried out by NASA's own scientific staff, and also, in collaboration with 
scientists in other institutions. We hope this description and evaluation of 
the model's performance will «nable the individual researchers to judge its 
utility as a GCM-experlmentation tool. A bibliography of works dealing with 
the present and the earlier versions of the GLAS climate model Is enclosed as 


an Appendix 


2. Model description 


The GLAS cllmete model has evolved from the GISS GCM, described by Somaer- 
vllle et al » (1974) and Stone et al . (1977), which in turn, was based on the 
3-level 1970 version of the GCM developed by A. Arakawa (1972) and Y. Mints at 
UCLA. The model uses a form of the o coordinate system (Phillips, 1957), with 
nine layers, all of the same a thickness, and a 4 x 5 degree latitude-longitude 
grid, which is modified in the polar regions as discussed below. The variables 
are staggered in the horixontal according to the B-grid scheme of Arakawa and 
Lamb (1977), and in the vertical according to the scheme of Lorenz (1960). The 
upper boundary of the model is at 10 mb. The prognostic atmospheric variables 
are the surface pressure, the zonal and meridional components of the horizontal 
wind, the temperature, and the water vapor mixing ratio. The prognostic boundary 
variables are the bulk ground temperature and ground wetness, and the snow depth. 

Although the current model is similar to the GISS GCM presented by Somer- 
ville et al. (1974) and Stone et al . (1977), significant changes have been made 
in both the finite-difference schemes and the physical parameterlzatlons. As 
described by Halem and Russell (1973), Halem et al. (1979), and Herman and 
Johnson (1978), the model incorporates a "split grid," in which the number of 
grid points on a latitude circle is systematically reduced near the poles. This 
allows the use of a ten-minute time step, with only weak longitudinal smoothing 
at high latitudes. The idea is similar in principle to that proposed by Kurihara 
(1965), but is applied at only a few interfaces separating latitudinal bands 
within which the horizontal resolution is uniform. The resolution in each lati- 
tudinal band is twice that of its poleward adjacent neighbor. The differencing 
scheme maintains the quasi-conservative properties of kinetic energy and mean- 
square vorticity, as in the Arakawa scheme. Although written to handle an 
arbitrary number of bands, the geometry used in the present model has five 
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latitudinal banda, conaiating of a central region extending from -66* to <f66*y 
bounded by one set of banda from 70* to 78*, and a aecond aet from 82* to 86*. 
The latitudinal reaolutlon la 4* in ali bands. The longitudinal resolutions 
are 5*, 10*, and 20*, respectively. 

Once every simulated half hour, a sixteenth-order Shapiro filter (Shapiro, 
1970) is applied, in the longitudinal direction, to the sea level pressure, and 
to the potential temperature and the wind components on the a surfaces. The 
filter was introduced as an ad hoc device to suppress the "checkerboard” noise 
in the prognostic fields, which arises from the inability of the B-grld to sim- 
ulate the geostrophic adjustment process at the two grid interval level (Arakawa 
and Lamb, 1977). However, a substantial price has been paid for this noise 
suppression. Filtering the potential temperature on the a surfaces leads to a 
systematic cooling of the air over mountains and a systematic warming in the 
neighboring valleys. Filtering of the winds dissipates kinetic energy. Filter- 
ing of the sea-level pressure Interferes with every conservation property of 
the model (except mass conservation). The magnitudes of some of these false, 
filter-induced sources and sinks have been determined. The filtering of poten- 
tial temperature generates 2.5 W m~^ of potential energy, and filtering of the 
winds dissipates 1.5 W m~^ of kinetic energy. 

The Matsuno forward-backward time differencing scheme la used. This scheme 
tends to damp high frequencies, and so It helps to control computational noise 
In the model, but at the cost of almost doubling the computing time required 
for the dynamics. 

We use the cumulus parameterization developed by Arakawa (1969; see also 
Haltlner, 1971) for the three-level UCLA GCM; but as modified for use In a 
nine-level GCM by Sommervllle et al. (1974; see also Helfand, 1979). Although 
this Arakawa cumulus parameterization predates the cumulus parameterization 
theory of Arakawa and Schubert (1974) and Lord and Arakawa (1981), t.any key 
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concepts eaployed by Arekewe and Schubert were already present In the 1969 
paper. Including the cuanilus mess flux, a spectrum of clouds distinguished by 
their fractional entrainment rates, and closure formulated in terms of the 
stabilisation of the environment by the clouds. Helfand (1979) modified the 
parameterisation to include very deep clouds and cumulus friction, but these 
fiodifications are not incorporated into the present model. 

i 

The model also includes latent heat release due to large-scale saturation, 
which occurs when the relative humidity exceeds lOOX. 

The short-wave radiation parameterization is that presented by Lads and 
Hansen (1974). It includes absorption by ozone, water vapor, and clouds. 

The long-wave radiation parameterization of the current model is based on 
the method of Wu (1976), Wu et al. (1978), and Wu (1980). The parameterization 
includes a water vapor transmittance chat uses a statistical band model with 
the strong line version of the Curcis-Godson (Godson, 1955) approximation; the 
water vapor dimer effect in the 8-13 window region; line-by-line precalcula- 
tion of CO 2 transmittance including fundaomntal bands, hot bands, and isotopes; 
tables of ozone transmittance calculated by Dr. N. Scott (private communication); 
a special treatment of the nearby layer quadrature for increased accuracy in 
the numerical integration; and the effects of clouds. In order to calculate 
the incoming long wave flux at the 10 mb level, the top of Che model, zonally- 
averaged climatological temperatures at 1 mb and 5 mb are prescribed. In order 
to save computing time, the time step for the long-wave radiation is 5 hours 
for the simulations presented in this paper. 

Clouds are assumed to occur if and only if the model predicts cumulus con- 
vection (restricted to the lowest six layers) or large-scale saturation (in any 
layer). No parameterization of subgrld-scale fractional cloudiness is attempted; 
both supersaturation and convective clouds are assumed to completely fill a 
grid box. Moeng and Randall (1982) have recently studied the sensitivity of the 
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Bodel to this sxtrsoa psrsMtsrissticn of convectlvs cloudlnsss* Ths optlcsl 
propsrtles of ths clouds sre proscribed ss shown in Tsbls 2.1* In ths tsrrss^ 
trisl rsdistion psrsMtsri*.stion» oil clouds srs sssuasd to bs block bodies. 

Ths plsnstsry boundsry Isysr (PBL) persasterisstion wss foraulstsd by 
Kstaysms, (ss reported by Arsksws, 1972) end aodifisd by Sosaerville, st al. 
(1974). We hev« revised the pereaeterisstion in order to correct e defect 
noted by Chamey et si. (1977): earlier versions of the aodel produced too 

little evaporation over land and too much over the oceans* We found a high- 
frequency (in time) coaput^tlonal node in the surface evaporation and sensible 
heat flux fields. A simple revision of the nunericsl algorithm for computing 
these fluxes eliminated the computational mode, and the revised aodel produces 
more realistic surface flux distributions (Sud and Ab; les, 1981). 

The model incorporates a prognostic soil aoisture, following the UCLA 
model (Arakawa, 1972). but with different forwilaticns for the runoff, snow 
melt, and potential evapo transpiration, as developed by Lin et al. (1978). 

A comprehensive, detailed, and up-to-date documentation of the aodel and 
its post-processing routines is in preparation as a NASA Technical Memorandum 
(Sud et al., 1982). 
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Table 2.1 Prescribed cloud properties in GLAS solar radiation calculations. 
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3. Boundary and initial conditions 

The present GLAS climate model uses the externally Imposed boundary con- 
ditions listed In Table 3.1. The observed climatological values of these 
boundary conditions are available only In the form of monthly means. A contin- 
uously changing set of prescribed boundary conditions was obtained by assigning 
the observed values to the midpoint (In time) of the appropriate month, and 
computing the values on all other days by linear interpolation. Table 3.1 
gives the source of the data for the winter and summer slnnilations, and the 
name of the GLAS data set containing the data. Fig. 3.1a shows the topography 
used for both runs, %rhile Figs. 3.2 and 3.3 show the sea surface temperature, 
land-sea and snow-ice distributions, and surface albedo for the months of 
February and July. 

As revealed by a comparison of Figs. 3.1a and b, the topography has been 
greatly smoothed, and the land-surface heights for Greenland substantially 
reduced. This was necessary in order to reduce the level of computational 
noise in the earlier versions of the model. 

Observed values of the surface albedo are prescribed for the land, ocean 
and desert regions are shorn in Figs. 3.3a and 3.3b. This data was prepared 
from that described by Posey and Clapp (1964). Although the snow depth is 
predicted by the model, prescribed albedoes were used exclusively in the present 
runs. The snow and ice distributions used to determine the surface albedo are 
those listed in Table 3.1 and plotted in Fig. 3.2. 

The winter simulation described in this paper started with the observed 
state of the atmosphere for OOZ on January 1, 1975. The summer simulation is 
started with the observed state of the atmosphere for OOZ on June 15, 1979. 

In both runs, the initial ground temperatures for land and sea ice points were 
assumed to be equal to the suface air temperature. The initial soil moistures 
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were obtained from climatological surface 
by Stone et al, (1977). The Initial snow 
climatological data of Matson (1978), and 
as functions of the surface albedo. 


air relative humidities, as explained 
boundaries were obtained from the 
the Initial snow depths were specified 
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Table 3.1 


Prescribed 
Boundary Condition 


Topography 

Sea Surface 
Temperatures 


Surface 

Albedoes 


Snow/ Ice 
Distributions 


Data 

Source 


NASA/ Goddard Institute 
for Space Studies 

NOAA/ Geophysical Fluid 
Dynamics Laboratory 


Posey and Clapp ( 1964) 


British Met. Office (1977) 
Matson (1978) 


GLAS 

Itota Set Name 

C1APF79T 

C1APF79T 

C1APF79A 

C1APF79S 
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Figure 3.1a Topography used in the winter and summer simulations. The con 
tour interval is 500 m. 
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Figure 3.1b Topography obtained by simple area averaging on to a 4 x 5 degree 
latitude-longitude grid. The contour Interval is 500 m. 



Figure 3.2a February sea surface temperatures, land-sea distribution, an 
snow/ ice distributions. The symbol meanings are: L = bare 

land, S = snow-covered land, 1 = ice. 





Figure 3.2b July sea surface temperatures, land-sea distribution, and snow/ 
ice distributions. The symbol meanings are; L * bare land, 

S » snow-covered land, 1 = ice. 
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4. Mean prognostic fields 
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a. Sea level pressure. Figures 4.1a and 4.2a show the 16 year mean sea 
level pressure for January and July reproduced from Godbole and Shukla (1981), 
and Figures 4.1b and 4.2b show the raoflel simulated mean sea level pressure for 
the single February and July, respectively. 

For February, the simulation of the prominent Northern Hemispheric circul- 
ation features, such as the Aleutian low, the Icelandic low and the Siberian 
high are fairly realistic. When comparing a single month model simulation with 
the mean climatology, one must keep in mind that there is considerable inter- 
annual variability in the observed monthly means. For example, the structure 
and Intensity of the Aleutian and the Icelandic lows show considerable varia- 
bility over the 16 years analyzed by Godbole and Shukla (1981) and the model 
simulation that is shown here falls well within the range of this variability. 
However, the discrepancy in the structure and intensity of the Siberian high 
(Compare Figs. 4.1a and 4.1b) is too large to be accounted for by the interannual 
variability or by differences between January and February. It has been conjec- 
tured by Halera et al. (1979) that the lack of ^ well organized Siberian high is 
due to excessive stratus formation over the region. Tne most serious deficiency 
of the simulated sea level pressure field in February is in the Southern Hemi- 
sphere. The bserved sea level pressure is zonally uniform between about 40°S- 
60°S, while the simulated pressure field shows a large eddy structure in this 
zone. It has been our general experience that the simulation of the Southern 
Hemisphere pressure field is very sensitive to the model physics and numerics* 

For July, the simulated sea level pressure field is also in fairly good 
agreement with the observations. The high pressure cell over the North Atlantic 
is very well simulated both with regard to intensity and to structure. Over the 
north Pacific the simulated high pressure cell is too small in the east-west 


direction, and Its major axis Is oriented NNW-SSE Instead of WSU-ENE. The 
Asiatic monsoon low Is well simulated. The subtropical high pressure cells In 
the Southern Hemisphere are well simulated over the Indian Ocean, Atlantic 
Ocean and the eastern Pacific but not over Australia. The westerly zoasl geo- 
strophlc flow that Is sho«m by the Isobars between 40*S and 60” S Is fairly well 
simulated over the south Atlantic and south Indian Ocean; but In the south 
Pacific Ocean the simulated geostrophlc flow Is easterly between 40*S-‘60”S. 

For summer and winter, the model shows a common deficiency of distinct non** 
zonal Isobar configurations ^r\ the southern hemispheric mld**latltudes. 

b. Geopotentlal heights. Figures 4.3a and 4.4a show the observed (Crutcher 
and Davis, 1969) 200 mb geopotentlal height field for winter (December, January, 
February) and summer (June, July, August), and Figures 4.3b and 4.4b show the 
corresponding model simulation results for February and July respectively. 

For February, the model simulation of the trough over the northeast USA 
and eastern Canada, the ridge over the eastern Atlantic and the Jet stream over 
Japan are quite realistic. The ridge over the west coast of the USA is displaced 
to the east. The simulation of the upper level antlcyclonic circulation over 
the tropics is reasonable. The actual values of the geopotentlal, however, are 
higher than the mean climatology. The southern hemispheric flow Is also well 
simulated. Model simulations show more distinct eddy structure in the southern 
hemispheric circumpolar flow along 40“-60"S. It will be shown in section 5 that 
the model simulated stationary variance is in excellent agreement with the 
observed stationary variance. 

For July, the model simulated 200 mb flow shows large amplitude short waves 
over North America and the adjacent Pacific Ocean. The ridge over central U.S.A. 
is well simulated but the trough off the west coast of U.S.A. is not present 
In the observed climatology. Although the local maximum associated with the 


Tibetan high la not wall alimtlated, the large scale anticyclone over the Asiatic 
nonsoon Is reasonable. The geopotential halght in the tropical belt Is higher 
by about 200 gpm than cllnatology. This Is due to a gradual warning of the 
tropical atmosphere In the nodel slnulatlon. We do not yet know the reason 
for this problem. The location of the simulated mid-Atlantic trough does not 
agree with the observations. The observed and simulated upper-level flows are 
both predominantly sonal in the Southern Hemisphere. 

The surface circulation during July and the upper level circulation during 
January are, In general, better simulated than the surface circulation during 
winter and the upper level circulation during summer. 

c. Flow at 200 mb. Figures 4.5a and 4.6a show the observed 200 mb flow 
pattern for February and July, and Figures 4.5b and 4.6b show the correspond- 
ing model simulations for February and July, respectively. The observed 200 mb 
flow for February was obtained from the NMC analysis of February, 1979 and for 
July it is reproduced from Sadler (1975). 

For February, the locations and the intensities of the Jet streams in both 
hemispheres are well simulated. The model simulated jet stream speeds in both 
hemispheres, are, in general, stronger than the observations of February, 1979, 
especially in the Southern He<nlsphere. Away from the equator, the flow patterns 
are quite consistent with the geopotential height fields shown earlier. 

The simulated July 200 mb geopotential height field does not compare well 
with the observations. One deficiency is the weakness of the summer monsoon 
easterly jet stream along 10“N. The subtropical jet stream near Japan is 
reasonable in terms of location and intensity, but the flow pattern over the 
subtropics of the Western Hemisphere is not realistic. 

d. Zonal wind. Figures 4.7a and 4.8a (solid lines) show the observed xonal 
wind tor February and July 1979, and Figures 4.7b and 4.8b (solid lines) show 
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the corresponding leodel simulations for February and July» respectively* 

For February, the most conspicuous deficiency of the simulation is the 
absence of the closed maximum near 200 mb which is seen in the observations* 

The model simulated zonal wind increases with height and tilts cowards the north 
so that the model simulated zonal wind near the upper boundary is stronger than 
is observed. This model deficiency la related to very low model temperatures 
in the the upper tropospheric polar regions. The locations of the strongest 
zonal wind maxima in both hemispheres Is fairly well simulated. The zonally 
averaged equatorial easterlies are a little too weak. The configuration of the 
zero wind line is not unreasonable except that it is too low south of the 
equator. 

The dashed lines in Figures 4.7a and 4.7b show the zonally averaged temper- 
ature. The largest discrepancy is found near t'ne upper boundary over the polar 
region';. Dynamical heat fluxes are not sufficient to compensate for the exces- 
sive radiative cooling associated with high level supersaturation clouds which 
are treated as black bodies. The middle troposphere over the polar regions is 
warmer than the observations. The model simulated tropical atmosphere is consi- 
derably warmer than the observations. In a separate study Shukla and Wallace 
have found that most of the warming of the tropics occurs in first 30 days. 

For July, the observed and simulated zonal wind fields are In very good 
agreement, except that the Jet in the Southern Hemisphere increases monotonlcally 
with height so that the maximum value occurs near the upper boundary. Zonally 
averaged temperatures (dashed line) near the northern pole are lower in the 
model simulations than In the observations. The model simulated tropical atnu- 
sphere is aga^n warmer than observations and this may be one of the reasons for 
a stronger westerly Jet in the Southern Hemisphere. 

e. Zonally averaged specific humidity. Figures 4.9a and 4.10a show the 

observed zonally averaged specific humidity for February and July, 1979, and 
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Figures 4.9b and 4.10b show the corresponding model simulations for February 
and July, respectively. For both seasons, the general latitudinal and vertical 
structure of the humidity field quite reasonable. However, the simulated 
specific humidity at the lower levels Is larger than observed. As will be 
shown later, the simulated rates of evaporation and precipitation are comparable 
to the observations. Therefore, the high values of specific humidity may be 
related to the boundary layer parameterization and the values of surface transfer 
coefficients used In the simulations. 
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Fi;?ure 4.1a Observed mean sea level pressure (-1000 mb) for January (Godbole 
and Shukla, 1980). 



Figure 4.1b Model simulated sea level pressure (-1000 mb) for February 
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Figure 4,2b Model simulated sea level pressure (-1000 mb) for July 
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Figure < 4 . 3b Model simulated geopotential height (gpm) at 200 mb for February. 








Figure 4.4a Observed mean geopotential height (gpm) at 200 mb for June, July 
and August (Crutcher and Davis, 1969). 
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Flgura 4.4b Model simulated geopotential height (gpm) at 200 mb for July 
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Figure 4.5a Observed wind direction and speed (ras“^) at 200 mb for February, 
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Figure 4.5b Model simulated wind direction and speed (ms ) at 200 mb for 
February. 
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Figure 4.7a Observed zonally averaged zonal wind (ras~^, solid lines) and 
temperature (“K, dashed lines) for February, 1979. 




Figure 4.7b Model simulated zonally averaged zonal wind (ms~^ , solid lines) 
and temperature (°K, dashed lines) for February. 














Figure 4.10a Observed zonally averaged specific humidity (g kg *) for July, 
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Figure 4,lub Model simulated zonally averaged specific humidity (g kg 
for July. 
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S* Fluctuations in space and time 

a. Energetica-hemiapheric energy cyclea* Hie energy cycle in an entity 
which summarizes the energy content of the atmosphere, its resolution into 
various components, and the transformations between these components (Lorenz, 
1955; Oort, 1964). Figures 5.1 and 5.2 present the energy cycles of the simu- 
lated winter and summer seasons, and those estimated for the real atmosphere 
(Oort and Peixoto, 1974). Each energy cycle is computed In the space-time 
domain (Oort, 1964) and on a hemispheric basis^. 

The observed energetics for January in the Northern Hemisphere (Figure 
5.1a) Indicate a strong (baroclinic) flow of energy from zonal available poten- 
tial (Pm) to eddy available potential (Pg) to eddy kinetic (Kg) energy, and a 
weaker (barotropic) flow from eddy kinetic to zonal kinetic (Km) energy. Also 
indicated is a conversion from Km to ?m, although this estimate has a large 
degree of uncertainty associated with it (Oort and Peixoto, 1974). These 
energy transfers are reproduced remarkably well in the model Northern Hemisphere 
energy cycle for February (Figure 5.1b). The sizable conversion of Km to Pm 
in the model suggests that the energy conversion of the thermally indirect 
Ferrel cell overcompensates that of the thermally direct Hadley cell. This is 
not inconsistent with the much larger mass flux associated with the Hadley cell, 
since the Ferrel cell exhibits a much larger temperature difference between the 
ascending and descending branches. The energy amounts in the model simulation 
agree well with the observations except for Km! the fact that the model Km is 
almost twice that of the observations is consistent with th^ model's unrealisti- 
cally strong subtropical Jet. Given the large uncertainties in estlmrtii.g the 

^ The model conversions from potential to kinetic energy were computed using 
the formulation, while the estimates given by Oort and Peixoto (1974) 

used the ”V*7Z'* formulation. Hence it was necessary to transform the observed 
estimates to the "(ua" formulation by making use of the observed fluxes of 
energy across the equator (Oort and Peixoto, 1974; Peixoto and Oort, 1974). 

BLANK 
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generation and dlsalpatlon of energy in the real atiaoephere (Oort and Peixoto, 
1974), it is difficult to assess the model's performance with regard to these 
quantities. The dissipation of kinetic energy in the model's Northern Hemisphere 
(in winter) is within a factor of two of the observed estimate, with the model 
showing less dissipation of eddy kinetic energy titan the observations. The 
model's winter generation of P^ (due to both dry and moist processes) is much 
stronger than the observed estimate. This may be due to excessive cooling near 
the winter pole and warming of the tropical atmosphere. The model also predicts 
a small net diabetic loss of Pg, in contrast with the estimate for the atmosphere. 

The Southern Hemisphere energy cycle for the winter (July) simulation is 
indicated in Figure 5.1c, and it Is seen to be quite similar to the simulated 
Northern Hemisphere winter cycle (Figure 5. lb). The only energy conversion which 
is substantially different is the barotropic conversion from Kg to K^, which is 
of opposite sign in the two hemispheres. The energy sources and sinks are also 
qualitatively similar, and the zonal available potential and eddy kinetic ener- 
gies are close in magnitude. Pg Is greater in the Northern Hemisphere, while 
Km is greater in the Southern Hemisphere. 

The observed summer Northern Hemisphere energy cycle (July) is presented 
in Figure 5.2a. The barocllnlc energy conversions, P^ to Pg and Pg to Kg, still 
have the same sense, although they are smaller than in January. The weak (baro- 
tropic) conversion from to Km Is again present, but the estimate of the con- 
version between Km and Pm has changed sign from the January observations. The 
amounts of the different types of energy are considerably less in summer than 
in winter. The energy conversions of the model Northern Hemisphere summer (July, 
Figure 5.2b) are In the same sense as observed, but they are generally of much 
larger magnitude. This is particularly true of the conversion from Pg to Kg. 
Furthermore, the amounts of energy in the model simulation are somewhat larger 
than observed, the discrepancy being the greatest for Pm. Frictional dissipation 
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In the Bodel is epproxinately the sasM as the obeerved eatlnate, hut the total 
model generation of Pm and Pg greatly exceeds the eatimatea of the real atmo- 
sphere. In general » the slwilated Northern Hemisphere summer circulation 
appears to be ov«rly vigorous. 

Comparing the model energy cycles for Northern Hemisphere suauaer (July, 
Figure 5.2b) and Southern Hemisphere sumer (February, Figure 5. 2c), It Is seen 
that the conversions are generally similar In magnitude and are In the same 
direction. It should be noted, however, that the conversions Pm to Pg and Kg 
to Km are larger In the Southern Hemisphere sususer. Both hemispheres contain 
similar amounts of Kg and Pg, but the Southern Hemisphere % and ?m are signi- 
ficantly greater than those of the Northern Hemisphere. The sources and sinks 
of energy are comparable In the two summer hemispheres, although the generation 
of Pg by dry dlabatlc processes Is of opposite sign. 

Meridional energy transports. This subsection presents the meridional 
transports of the various components of energy density, viz., sensible heat 
(Internal energy), potential energy, kinetic energy, and the latent energy of 
condensation of water. These transports are time-averaged and vertically inte- 
grated, and are also divided Into eddy and zonal-mean components. 

The observed zonal-mean transports of the various forms of energy density 
are presented for February In Figure S.3a, while the mean model transports are 
given In Figure 3.3b. The source of all the observed transports reported here 
Is Oort and Rasmusson (1971). The observed mean transport of potential energy 
Is directed opposite to the transports of latent and sensible heat; potential 
energy Is carried northward by the Hadley cell, whereas sensible and latent 
heat are transported southward. North of 26N, the presence of the Ferrel cell 
causes this pattern to reverse, although the mean latent heat transport becomes 
very small. Because the transports of the Individual components tend to oppose 


•ach other, the total energy tranaport la relatively avail. It la aleo con- 
aiatently In the direction of the potential energy tranaport. TIm vean trana- 
port of kinetic energy la not depicted, elnce It le too evall to play any role 
In the overall energy tranaport. The February zonal-flwan energy tranaporta of 
the node! are given In Figure 5.3b. The baeic pattern of tranaporta agreea 
well with the obeervatlone, with the potential energy tranaport again oppoaing 
that of sensible and latent heat. However, the transition froa northward to 
southward flow of total (or potential) energy occurs about 8 degrees too far 
north. In addition, the magnitude of the sensible heat and potential energy 
transported by the model's Hadley cell is too small, as Is the amount of sen- 
sible heat flow In northern mid-latitudes. (Note the difference In scale 
between Figures 5.3a and 5.3b.) The Hadley cell transport of latent heat agrees 
well with observations. 

The observed eddies transport both latent and sensible heat northward in 
February, the total northward energy flow reaching a broad maximum In northern 
mid-latitudes (Figure 5.3c). The eddy transports of potential and kinetic energy 
are both quite small. Although the model results for the latent and sensible 
heat carried by the eddies (Figure 5.3d) show basic agreement with the observa- 
tions north of 30N, there are discrepancies. The magnitude of the total model 
transport Is somewhat too large, and the peak In sensible heat flow occurs 4 
degrees (one grid point) too far south. Between 30N and the equator, the model's 
eddies transport sensible heat southwards, while the latent heat transport 
remains near Its peak value. Neither of these features Is observed. In the 
Southern Hemisphere, the (poleward) eddy transports of latent and sensible heat 
resemble those In the northen hemisphere, although the amount of sensible heat 
carried Is considerably smaller. The potential and kinetic energy transports 
are again quite small. 
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It li quite epperent fron Figures StSe-S.Sd thet elaoet ell of the poten- 
tiel energy treneport ie eccoapUehed by the aeen flow. In contreet, Figure 
9.3e indicates that in the cese of kinetic energy, it ie the eddiee tdiich doai** 
nets the trensport. The aodel’s eddy kinetic energy flow egrees with the 
observed only in the grossest sense: the trensports ate both to the north froa 
about 14N to 62N. Not only does the aodel trensport auch aore energy then is 
observed, but It trensports kinetic energy southwerd north of 62N, in contrest 
to the observations. Also, the double aaxlBua in the observed transport is 
absent in the model results. 

The observed zonal-aean transports of energy in July (Figure 5.4a) strongly 
resemble the winter transports, with the Hadley cell now shifted southwards end 
directed towards the Southern Hemisphere. The mean energy transports of the 
model in July (Figure 5.4b) Indicate the sane deficiencies present in February, 
namely that tho transports of sensible heat and potential energy are too smell. 
(Note the difference in scale between Figures 5.4a and 5.4b.) Also, the latitude 
of peak transport of the Hadley cell is located at least 4 degrees too far 
north. Gauging the strength of the model winter hemisphere Ferrel cell by its 
transport of potential energy, one concludes from Figures 5.3b and 5.4b that 
the cell is weaker in July than it is in February. 

The observed July eddy transports (Figure 5.4c) are naturally much weaker 
than in February, with the peak magnitude of northward latent heat transport 
greater than that of the sensible heat transport. This overall pattern is 
reproduced in the Northern Hemisphere July model results (Figure 5.4d), although 
the model peak of northward latent heat transport is too large by almost a 
factor of two. The position of the peak In sensible heat transport is slightly 
too far south in the model, although the peak magnitude is correct. Also, the 
transition from northward to southward energy transport occurs near the equator, 
as is observed. The flows of latent and sensible heat in the model Southern 
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(wlnt«r) IteHisptmre gen«r«lly reieabl* those of the February northern (winter) 
Heaisphere (Figure 5* 3d), with the exception that the July sensible heet trans- 
port is a factor of two leas than the February transport. 

The July eddy transports of kinetic energy, both as simulated and observed, 
are shown in Figure 3.4e. The only significant observed feature in the Northern 
(sumser) Heaisphere is a peak of northward transport at 34N. This is absent in 
the model results, which indicate weaker maxima of northward transport at 42N 
and 58N and maxima of southward transport at 28N and 66N. In the SMithern Hemi- 
sphere, the model's eddies carry a great deal of kinetic energy southward near 
32S, and, further south, an equally large amount northward. Whether the very 
sizeable convergence of eddy kinetic energy at 40S is realistic or not Is a 
question which cannot be answered with the data presently in hand. 

b. Stationary waves in the geopotential height field. The time-averaged 
departures from zonal symmetry of the basic meteorological fields reflect the 
lack of such symmetry at the lower boundary. That is, the presence of station- 
ary waves is due to the existence of orography and geographically fixed, asym- 
metric, heat sources. Thus, in order for a GCM to reproduce the observed 
three-dimensional structure of the stationary waves, the model must simulate the 
geographically fixed heat sources and flow over orography in a realistic manner. 

In this section results are presented for the stationary (time-averaged) 
component of the geopotential height, from which the stationary temperature and 
wind fields may be estimated by using hydrostatic and geostrophic balance respec- 
tively. Figure 5.5a shows the observed winter variance around a latitude circle 
of the stationary geopotential height, summed over the planetary waves (wave- 
numbers 1-4). The term "stationary” refers to a 90-day average commencing on 
January 1. This average was computed from NMC analyses for each of the years 
1963-1977, and the mean over these years was taken. The correaponding stationary 
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planetary wave variance for the winter GCM ainulation la preaented in Figure 
5.5bl. ’’Stationary" here refera to a tlae-average over the flrat 90 dayit of 
the aiiiulation, which waa initlallted from data valid for January 1, 197S. The 
■odel reaulta ahow good overall agreeiMnt with the obeervatione, with both 
variances reaching a naxiaun in the aid-latitude upper troposphere. The model 
variance is slightly too weak up to the 250 ab level, above which it fails to 
ahow the decrease with height apparent in the observations. This failure may 
be related to the unrealistically low polar teaperatures and strong zonal winds 
at the upper levels. The model stationary planetary wave variance in the 
Southern (suamer) Hemisphere (Figure 5.5b2) is almost an order of magnitude 
smaller than in the Northern Hemisphere, due to the relative lack of orography 
and geographically fixed heat sources in the Southern Hemisphere. It should 
be noted that the zonal variance of the stationary height field in the lowest 
four wevenumbers is (in winter) several orders of oMgnicude larger than the 
variance in higher wavenumbers, indicating that the total stationary variance 
will strongly resemble the planetary wave variance. 

A more detailed description of the stationary waves Is afforded by a zonal 
wavenumber decomposition of the geopotential height field. Figures (5.6a)- 
(5.6f) present the amplitude and phase^*^ of wavenumbers 1-3 (separately) for 
both the observed and simulated winters. Comparison of Figures (5.6a) (15 year 
mean observations) and (5.6b) (CCM) shows that the amplitude of wavenumber 1 
in the GCM agrees well with the observations in that both amplitudeji reach a 


2 The phase is defined as the relative position of the ridge, i.e. the phase 
is defined by writing the wave as A cos (kX - ^), where k is the zonal wave- 
mmiber and X the longitude. 

3 The plots of amplitude and phase were constructed from data at only five vertical 
levels (200 mb, 300 mb, 500 mb, 700 mb and 850 mb). This vertical resolution was 
not adequate to resolve certain ambiguities in the behavior of the phase as a 
function of latitute and height. Certain portions of the phase plots hence were 
left blank. 
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maximum in the mid-*latitude upper troposphere. However, the GCM amplitude does 
not decrease above 300 mb as does the observed amplitude, and the maximum in 
the GCM amplitude is both somewhat too strong and located about 6 degrees too 
far north. The phase of the GCM's wavenumber 1 tilts westward with height in 
mid-latitudes, but with less slope than is observed. The eastward phase tilt 
with Increasing latitude evident in the observations is almost completely 
absent in the model. Further, the actual values of the GCM phase are about 30 
degrees less than (to the west of) the observed values near the mid-latitude 
amplitude maximum. 

The behavior of the model wavenumber 2 (Figures 5.6c-5.6d) presents a con- 
trast to the wavenumber 1 results in that the phase shows better agreement with 
the observations than does the amplitude. The simulated amplitude is only half 
as large as the observed throughout mid-latitudes, although the overall latitude 
height structure is similar, with a broad mid-latitude upper tropospheric maxi- 
mum and one in the sub-tropical upper troposphere. Both the simulated and 
observed wave phase show a westward tilt with height and with (increasing) 
latitude in raid-latitudes, although the model phase tends to tilt too rapidly. 

In the vicinity of the raid-latitude amplitude maximum, the phase of the GCM 
wave is about 30 degrees too large (too far to the east). 

The amplitude and phase of the observed and simulated wavenumber 3 are pre- 
sented in Figures (5.6e)-(5.6f). The model amplitude agrees well with the 
observations, both with regard to the overall latitude-height structure and to 
the maximum amplitude, although the GCM maximum is located about 4 degrees too 
far north. In the vicinity of the mid-latitude maximum both the GCM phase and 
the observed phase show a moderate westward tilt with height and eastward tilt 
with (increasing) latitude. The values of the phase are also in rough agreement 
To the north of 60N, the observed phase tilts westward with increasing latitude, 
while the model phase continues to shift eastward with latitude in the lower 
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troposphere. In the subtropics • the nodel captures the observed westward phase 
shift with height but fails to capture the observed eastward shift with latitude 
in the upper troposphere. 

The observed suaner variance of the stationary coaponent of the geopoten- 
tial height, suaaed over wavenuabers 1-4, is presented in Figure 5.7a. The cor- 
responding GGM variance is shown in Figure S.7bl. These quantities were computed 
in the same manner as for the winter season, except that the averaging period 
was defined to be the 90 days starting on June 15. (This tiae period was 
chosen to correspond to the summer simulation, which was initialized with data 
valid for June 15, 1979.) The simulated planetary wave variance agrees well 
with the observations with regard to the location and strength of the three 
maxima: one in mid-latitudes at about 300 mb, and two in the subtropics, at 

830 mb and 200 mb. Discrepancies include the fact that the simulated mid- 
latitude variance does not decrease above 300 mb, the slight northward displace- 
ment in the GCM of the subtropical loaxiroum at 850 mb, and the weakness of the 
simulated subtropical maxlraura at 200 mb. It is noteworthy that the Southern 
(winter) Hemisphere .'.Itra-long wave variance (Figure 5.7b2) is now of the same 
order of !sagnltude as the Northern (summer) Hemisphere variance. 

The amplitude and phase of the stationary zonal wavenumber 1 in summer is 
shown in Figures 5.8a for the observations and in Figures 5.8b for the GCM. 

The latitude-height structure of the observed amplitude is very similar to that 
of the total planetary wave variance; one mid-latitude maximum in the upper 
troposphere is present, as are sub-tropical maxima at 850 mb and 200 mb. The 
model amplitude shows good agreement with the observations, the major discrep- 
ancies being the excess strength and southward displacement in the model of 
the midlatitude maximum. The GCM phase strongly resembles the observed south 
of 34N, exhibiting strong eastward tilt with height but little tilt with lati- 
tude. In mid-latitudes, however, the model phase is less realistic. The 
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observations show a slight eastward tilt with height and a stronger eastward 
tilt with (increasing) latitude south of the mid-*latltude oaxiauuB, and a moderate 
westward tilt with height and no tilt with latitude Just north of the maximum. 

The phase in the GCM, in contrast, tilts westward with height and eastward with 
(increasing) latitude through much of the raid-latitude region. 

The observed summer amplitude of stationary wavenumber 2 shows a good deal 
of structure (Figure 5.8cl); there are three distinct amplitude maxima and two 
distinct minima. The maxima again occur in the upper troposphere in mid- 
latitudes, and in the subtropics at 850 mb and 200 mb. The structure of the 
GCM amplitude (Figure 5.8dl) is generally similar, but differences of detail 
exist. The locations of the maxima are too far north, both upper tropospheric 
maxima are too strong, and the minima r.-e not as distinct as in the observa- 
tions. The behavior of the observed phase shows several regimes, and these 
are present to some extent in the model results. In the upper tropical tropo- 
sphere, the observed phase tilts little with latitude, and tilts westward with 
height. The model phase In this region tends to show much less shift with 
height and much more (westward) tilt with latitude. The observed phase in mid- 
latitudes tilts eastward with (increasing) latitude in the lower troposphere, 
and westward in the upper troposphere, up to the latitude of the amplitude 
maximum. The phase tilts eastward with height throughout this region. North 
of the amplitude maximum, the phase slopes westward with latitude, but shows 
little shift with height. The phase of the GCM wave reproduces this structure 
very well in an overall sense, but the values of the phase tend to be uniformly 
smaller (wave located too far to the west). 

The structure of the observed summer wavenumber 3 amplitude is quite simi- 
lar to that of the lower wavenumbers (Figures 5.8e), with three distinct maxima 
again present. In this case, however, the 850 mb maximum is located in mid- 
latitudes, at 38N. The GCM wave 3 amplitude again follows the pattern quite 
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well (Figures 5.8f), although the naxloa are located slightly too far north, and 
the upper-tropospheric mid-latitude maximum Is somewhat too weak. The observed 
phase slopes rapidly eastward with latitude In the subtropics, and slopes slowly 
westward with height. The model phase, however, shows a more rapid westward 
tilt with height and a westward tilt with latitude. Between 38N and 46N, the 
observations indicate almost no phase tilt with height and a strong eastward 
tilt with latitude, while the phase of the GCM shows a strong eastward tilt 
with height and no tilt with latitude. Between 46N and the mid-latitude ampli- 
tude maximum, both observed and simulated phases slope rapidly eastward with 
latitude, and very little with height. The region near the amplitude maximum 
is one with little phase change for either the GCM or the real atmosphere. 

c. Space-time structure of transient fluctuations-local analysis of variance. 
A useful description of the fluctuations in any basic field is afforded by the 
analysis of local variability discussed by Blackmon (1976). This method portrays 
the spatially local fluctuations of a particular field by the construction of 
maps of its RMS deviation in time. The variability due to low frequency motions 
or that due to higher frequency baroclinic activity can be separated by first 
filtering the data with a low-pass filter (retaining periods of 10-90 days), or 
a band-pass filter (retaining fluctuations of 2.5 to 6 days). Ir^ the applica- 
tion of this method, an estimate of the annual cycle is removed from the data 
before either the filtering or the computation of the RMS deviation. 

The RMS deviation of the (unfiltered) 500 mb geopotential height in the 
Northern Hemisphere winter is presented for the observations^ in Figure 5.9a 


^ Figures 5.9a, 5.10a and 5.11a are taken from Blackmon (1976). In that paper, 
the winter season was defined to be the period of 120 days commencing on Novem- 
ber 15. The data consisted of the NNC analyses, twice daily, for the years 
1963 through 1972, from which nine winter seasons were abstracted. The first 
four harmonics of the annual cycle were removed. The RMS deviations of the 
filtered and unfiltered fields were computed for each winter, and then averaged. 


- 55 - 


and for the GCM^ in Figure 5.9b. The observed aap shows a trend of increasing 
variability with higher latitude, at least up to mid-latitudes. Particularly 
noticeable are the regions of high variability in the northeast Pacific, in the 
north Atlantic, and over north-central Asia. The GCM results show the same 
general trend, although the values of the RMS deviation tend to be less. The 
observed high in the north Pacific is present in the model, although in much 
%raakened form, while the observed strong north Atlantic high is replaced by a 
pattern of weaker highs over Newfoundland and off the east coast of Greenland 
and a low over the northeast Atlantic. In addition, the model shows a region 
of relatively high variability over the west coast of the United States, a 
feature that is not observed. The level of noise in the GCM map relative to 
the observations Is to be expected, for the latter represents an average of 
nine individual atmospheric realizations, whereas the former represents a single 
(simulated) realization. 

The discrepancies between the simulated and observed RMS maps are further 
elucidated by the separate consideration of the low-pass and band-pass RMS 
maps. The low-pass filtered geopotential height field gets strong contribu- 
tions from blocking ridges, or positive anomalies that persist longer than, 
say 10 days (Blackmon, 1976). Their occurrence in nature is highly variable 
from one year to the next, so that the RMS low-pass field for any one year (or 
any one model simulation) may be quite different than the climatological mean. 
The average RMS low-pass map for the observations (Figure 5.10a) and for the 
GCM (Figure 5.10b) both strongly resemble their unfiltered counterparts. Thus 
the discrepancies noted between Figures 5.9a and 5.9b are in large part the 


^ Here the winter season is defined as the 90-day period commencing on January 1. 
The annual cycle was estimated by fitting each 90-day time series with a para- 
bola. The filters used are identical to those of Blackmon (1976). 
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differences between the clitsatology of persistent anomeltes and their presence 
or absence In one model simulation* 

The band-pass RMS, on the other hand, Is closely associated with regions 
having a high frequency of cyclonic activity, and so Indicates the location of 
the “storm tracks” (Blackmon, et al, 1977). A comparison of the observed band 
pass RMS (Figure 5.11a) with that of the GCM (Figure 5.11b) shows excellent 
agreement in terms of both the location and strength of the major areas of 
cyclonic activity in the north central Pacific and western Atlantic* 

The observed total RMS of the summer 500 mb height fleld^ Is presented In 
Figure 5.12a. The two oceanic regions of high variability In the winter map 
(Figure 5.9a) are again present, although they are weaker and their positions 
have been somewhat shifted. The simulated total RMS summer 500 mb height field 
(Figure 5.12b) is quite noisy compared to the observations. In the Atlantic, 
there is an extended region of high RMS, but the RMS is much too small* The 
GCM shows a rapid succession of highs and lows in the Pacific region where the 
observations indicate a broad high. The average RMS in this area is also much 
too small. 

The observed summer band-pass RMS height field (Figure 5.13a) shows the 
two oceanic maxima present on the winter map. The Atlantic maximum has been 
shifted to the east, and the magnitude of the summer RMS is generally much 
smaller. The corresponding (Northern Hemispheric) map for the GCM (Figure 5.13b) 
indicates that the Atlantic maximum has been realistically simulated, both with 
respect to position and to magnitude. The Pacific maximum, however, is too 
weak and is located too far to the west. 


^ The observed fields (Figures 5.12a, 5.13a) are taken from Blackmon (1976). 
The summer season was defined there to be the 123-day period from May 15 to 
September 14. The data base was identical to that mentioned in footnote 4, 
as was the methodology used to obtain the maps. 
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The RMS deviation for the 850 mb temperature field for the observations 
and for the model are shown in Figures 5.14a and 5»l4b, respectively^. The 
observations show greater variability over land than over oceans » especially at 
high latitudes, with no clear indication of the oceanic storm tracks. The most 
noticeable region of high variability is that over northwestern Canada, reflect- 
ing the alternation of high marine and low Arctic air temperatures (Blackmon 
et al. , 1977), Other (less striking) regions of high RMS are Greenland and an 
area in north-central Asia. The GCM map also shows more variability over land 
than over water. The maximum in the RMS over northwestern Canada and Greenland 
are not only present but are In fact unrealistically strong. In addition, the 
Atlantic storm tracks stand out more in the GCM than in the observations. The 
observed high over north-central Asia appears over eastern Siberia in the GCM, 
and also appears too strong. Both the observed and simulated low-pass KMS maps 
(Figures 5.15a and 5.15b) resemble the total RMS maps, with the exception that 
the model's Atlantic storm tracks now are not apparent. 

The observed band-pass RMS (Figure 5.16a) shows both the presence of high 
variability over the oceanic storm tracks and in the lee of the Rocky mountains, 
presumably reflecting the cyclogenetic activity there. The GCM band-pass RMS 
(Figure 5.16b) shows the presence of both the Pacific and Atlantic storm tracks. 
The lee-side Rocky mountain maximum does not appear on the simulated map. 

Wavenumber-frequency spectra. Space-time spectral analysis is an 
approach in which the longitude- and time-variance of a field is separated into 
components of distinct spatial scales (described by wavenumber) and time scales 
(measured by frequency). The calculations of the wavenumber-frequency spectra 

^ The observed RMS deviations are from Blackmon et al. (1977). The annual 
cycle was estimated from the 850 mb temperature field in that work by fitting 
a parabola to the time series at each grid point. Otherwise, the methodology 
is as described in footnote A. 
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reported here are described In Straus and Shukla (1981), to which the reader 
should refer for complete details^. 

The observed winter wavenumber^freciuency spectrum of the 500 mb geopoten- 
tlal height field at SON is shown in Figure S.17a. and the model results are 
given in Figure 5.17b. These plots s)k)w the spectral density as a function of 
zonal wavenumber and frequency (in cycles per day). Goth the observed and 
model spectra possess two noteworthy features: (i) the lowest frequencies and 

wavenumbers have the greatest spectral density, and (ii) a ridge of high spec* 
tral density extends diagonally from lower frequencies and wavenumbers to 
higher frequencies and wavenumbers. These two features are common to many of 
the spectra that are presented in this paper. The orientation of this ridge 
is very similar in the model and observed spectra, and is consistent (in a very 
approximate way) with the Rossby phase speed formula (Pratt, 1977). Hence, 
this ridge is referred to as the “propagation band" (Pratt, 1975). The major 
difference between Figures 5.17a and 5.17b is the relatively small magnitude 
of the model spectral density at low frequencies and low wavenumbers. This 
shortcoming of the model will be evident in further results (to be presented), 
and is also a falling of other GCM's (Pratt, 1979). 

The dashed lines in Figure 5.17a divide the wavenumber-frequency plane into 
four domains: (i) low frequency planetary waves (LFPW's), consisting of wave- 

numbers 1-A with periods of 7.5 to 90 days, (ii) low frequency synoptic-scale 
waves (LFSW's), consisting of wavenumber 5-10 with periods of 7.5 to 90 days, 
(ill) medium frequency planetary waves (MFPW*s), consisting of wavenumbers 
1-4 with periods of 2.8 to 6.9 days, and (iv) medium frequency synoptic-scale 

^ A brief description of the methodology is as follows: The winter (summer) 

season is defined as the 90-day period commencing on January 1 (June 15). 

The observational data consist of 15 years (1963-1977) of NMC analyses. Each 
observed season and model simulation had a separate parabolic estimate of the 
annual cycle removed from each grid point before the spectra were computed. 

The reported observational spectra have been averaged over 15 seasons. 
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waves (MFSW's), consisting of wavenumbers S~10 with periods of 2.8 to 6.9 days. 
The latter category Is associated with propagating, barocllnically active 
systems, while the LFSW category may be related to the growth and decay (life 
cycle) of such systems. There are numerous physical mechanisms capable of 
explaining the presence of the LFPW's in the atmosphere. Including the insta- 
bility of either axlsymmetrlc or non-axlsymiMtric mean states, the non-linear 
transfer of energy from synoptic-scale waves, or the propagation of essentially 
free modes. Which of these explanations Is valid Is not known. These categories 
will be used to compile the space-time spectral summaries presented in a later 
subsection. 

The observed spectrum of 500 mb geopotential height at 30N in winter is 
given in Figure 5.18a, and the GCM spectrum in Figure 5.18b. The observed 
spectrum at 30N is less Intense than at 50N, although the main features of the 
two spectra resemble each other. The simulated spectrum at 30N has the same 
overall shape as the observed, with the magnitudes of the spectral density being 
generally comparable. The model LFPW's are only moderately too weak, although 
the discrepancy is more striking for wavenumber 1. The simulated spectrum 
also shows a peak at wavenumber 6 and frequencies near 0.10 cpd, a feature that 
Is not observed. 

The observed (simulated) spectrum of 500 mb height at 50N In summer Is 
given In Figure 5.19a (5.19b), while the spectrum at 30N Is shown in Figure 
5.20a (5.20b). The observed summer spectra are much weaker than their winter 
counterparts, particularly at 30N. At 50N, the propagation band has been 
shifted to higher wavenumbers (consistent with the observations of Pratt, 1975), 
while at 30N it has ceased to exist altogether. The propagation band of the 
simulated summer spectrum at 50N does not appear to have shifted from its winter 
position, although It Is not very well defined. In addition, the shape of the 
model spectrum is quite different than observed, with the model low frequency 
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peak at wavenumbers 1 and 2 contrasting with the observed peak at wavenumbers 
A and 5. The model spectral density Is generally too small. At 30Nt the model 
spectrum does not show any signs of a propagation band, a feature which agrees 
wlnh the observations. But rather than Indicating one low frequency peak at 
wavenumber I, the model results show a series of peaks at wavenumbers 1, 3 and 
7. Here the low frequency model spectral density is much larger than Is observed. 

A further useful diagnostic defined In the space-time spectral domain Is 
the net propagation tendency, which measures the extent to which motions of a 
particular space and time scale propagate preferentially westward or eastward. 
(This measure is defined so that -1.0 corresponds to pure eastward propagation, 
■*■1.0 to pure westward propagation.) The net propagation of the observed (simu- 
lated) winter 500 mb geopotential height at SON Is presented in Figure 5.21a 
(5.21b), and that of the temperature® in Figure 5.22a (5.22b). 

Both the observed temperature and heights have a net propagation tendency 
of close to -1.0 In the MFSW region, consistent with the eastward propagation 
of synoptic-scale systems. On the other hand, it is clear that the observed 
low frequency planetary waves have little tendency to propagate preferentially 
in any direction; the values reported for those scales in Figures 5.21a and 
5.22a are small and positive for the height field, and small and negative for 
the temperature field. The observational study of Pratt and Wallace (1976) 
indicates that this may be due to the coexistence of two dominant low fre- 
quency planetary wave modes, one which propagates westward and one which propa- 
gates eastward. Although the westward propagating mode is dominated by the 
height field at 500 mb and the eastward propagating mode is dominated by the 

® The net propagation tendency of the temperature was taken from Pratt and 
Wallace (1976). The major differences between their methodology and that 
used in this paper are the definition of the winter season (they define it 
to be the 132-day period starting November 25), the choice of band width 
(they use 0.133 cycles per day as opposed to 0.0667 cycles per day used here), 
and the data base (they used only four years of NMC analyses). 


- 61 - 





tenpAtature field • both variables contribute to both nodeai leading to a certain 
amount of cancellation in the net propagation tendency. The observed pattern 
of net propagation tendency ia reproduced in the model results for temperature 
and geopotential (Figures 5.21b and 5.22b). Both fields show almost exclusively 
eastward propagation in the MFSW scales^ while the height (temperature) field 
shows a slight tendency to propagate westward (eastward) in the LFPW scale. 

The region of positive net propagation tendency in the model spectrum is. how- 
ever, not as extensive as in the observations. 

In order to provide a background for the space-time spectral analyses of 
heat and momentum fluxes to be discussed in subsections d and e. spectra of the 
u and v winds are presented here. The spectrum of the observed 300 mb u-wind 
(v-wind) in winter at 35N is presented in Figure 5.23a (5.24a), and the simu- 
lated spectrum (at 34N) in Figure 5.23b (5.24b). The observed^O u-wind spec- 
trum is similar to the spectra of the geopotential in that most of the spectral 
density is concentrated in the LFPW's, and in the appearance of a distinct 
propagation band. The propagation band of the model spectrum is somewhat dis- 
organized, although the LFPW maximum is well simulated. The spectral density 
is in general too low, especially for low frequencies at wavenumber 1. 

In contrast, the v-wind spectrum is very well simulated by the model. The 
propagation band is correctly represented as running from wavenumber 5 and lower 
frequencies to higher wavenumbers and frequencies. (That the low frequency end 
of the propagation band occurs at wavenumber 5 rather than 1 or 2 is a direct 
consequence of geostrophlc balance.) Also, the magnitude of the simulated 
spectral density is generally quite realistic, including the values at the max- 
ima. Howevever, the model spectrum does show a peak near wavenumber 7 and a 

The observed wind spectra are reproduced from Pratt (1977). The methodology 

of that paper is similar to that described in footnote 9. 
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frequency of 0*173 cpd, a feature that la not obaerved. (It la important to 
note that aome featurea may appear more localized in frequency in the GCM apectra 
than in Pratt 'a obaerved apectra because of the difference in band width alluded 
to in Footnote 9.) 

Space-time apectral aummariea. In order to obtain a repreaentation 
of the vertical and meridional structure of the wavenumber-frequency apectra, 
it la useful to be able to summarize each spectrum by a few parameters. It is 
for this purpose that the four space-time categories of the previous subsection 
(low frequency planetary waves, low frequency synoptic-scale waves, medium fre- 
quency planetary waves, and medium frequency synoptic-scale waves) were introduced 
The parameters to be discussed in this subsection are derived by integrating the 
spectral density over these four space-time regimes. 

The latitude-height structure of the observed winter LFPW variance is shown 
in Figure 5.25a. This variance reaches a maximum slightly above the 300 mb 
level at about 64N, and shows a decrease above this level. The GCM variance 
(Figure 5.25b) has a similar structure, except that the variance increases up 
to at least 200 mb. In addition, the model variance is somewhat too small 
throughout the northern hemispheric extratropical domain. 

The observed winter LFSW variance (Figure 5.26a) reaches a fairly sharp 
maximum at 300 mb near 46N. The structure of the simulated variance (Figure 
5.26b) is again similar to this, although the Jiaximum is much more diffuse in 
both the vertical and meridional directions. As before, the model variance is 
not large enough south of 62N, whereas north of this latitude it is too large. 

Figure (5.27a) presents the observed winter MFPW var:fance. Its structure 
is seen to be quite similar tc the LFPW variance, except that the maximum is 
pushed further northwards. The magnitude of the MFPW variance, on the other 
hand, is considerably less than that of the LFPW's. The deficiencies of the 
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ttodel MFPW variance (Figure S*27b) are quite aiailar to the deficlenciea of the 
model LFPW variance; the taodel variance la generally too amell* and doea not 
show the obaerved decrease above 300 mb. 

In contrast to these results, the sioulated NFSW variance (Figure S.28b; 
agrees quite well with the observations (Figure 5*28a), Not only is the model 
variance almost as large as is observed, but it also shows the observed decrease 
above the 300 mb raaxiraum. 

The observed summer variances in the LFPW, LFSW, MFPW and MFSW categories 
are presented in Figures 5.29a, 5.30a, 5.31a and 5.32a, respectively. In gen- 
eral, the latitude-height structures of the variance are quite similar to the 
corresponding winter results, except that the 300 mb maxima have moved somewhat 
farther north, and the magnitude of the summer variances is somewhat less. The 
summer LFSW variance does not fit this pattern, however, since the winter raaxl- 
mura at 300 mb has moved up to at least 200 tub in summer, and the magnitude of 
the variance has decreased only slightly. 

Figures 5.29b, 5.30b 5.31b and 5.32b give the GCM results for the variances 
In the LFPW, LFSW, MFPW and MFSW categories. All the model variances increase 
with height up to at least 200 mb, so that except for the LFSW's, their struc- 
ture is not in agreement with the observations. In addition, the model medium 
frequency (MFPW, MFSW) variances show indications of two maxima at 200 mb. The 
magnitudes of the planetary wave (LFPW, MFPW) variances are too small over the 
whole northern hemispheric extratropical domain, while In the case of the 
synoptic-scale wave (LFSW, MFSW) variances, the simulated magnitudes are too 
small in the vicinity of the observed maximum but are too large far to the 
south and north of this region. 

d. Northward eddy transport of heat. The analysis of meridional energy trans- 
sports presented in subsection a indicated that the total vertically integrated 


northward transport of heat by the eddies was fairly realistic In both thw 
winter and suaaer sli»ilatlons* In this subsectloni Che northward flux of heat 
Is diagnosed In auch more detail. The first part discusses latitude -height 
cross sections of both stationary and crnnslent tonally-averaged nerldlonsl heat 
fluxes » while the geographical structure of Che local transient heat flux and 
Its breakdown Into low-pass and band-pass cooponents are studied In the second 
part. Some representative space-time spectra are presented In the third part. 

Stationary and transient eddy heat flux. The observed (simulated) 
time mean February heat flux due to Che stationary eddies Is shown as a function 
of latitude and pressure level^^*^^ In Figure 5.33a (5. 33b). The agreement 
between the model results and "he observation Is quite good. The observed mid- 
latitude maximum of positive heat flux Is vrell reproduced, as Is the secondary 
maximum at 200 mb. Discrepancies Include the presence of a sizeable negative 
heat flux in the model subtropics at 200 mb (which Is not observed), the Incor- 
rect sign of the model heat flux north of 70N In the lower troposphere, and the 
slight equatorward displacement of the 200 mb maximum In the model results. 

The time mean heat flux due to the translert eddies in February Is given 
for the model^^ In Figure 5. 34b and for the observations in Figure 5.34a. It 
is important to note that no prefiltering of the data has been dune before com- 
puting the transient component, eo that this component includes a contribution 
due to the annual cycle. Both the observed broad ma 'imum of positive flux in 
the mid-latitude lower troposphere and the weaker maximum near 200 mb are well 
reproduced by the model. Several differences exist, however. The region of 

The observed results for the stationary and transient heat flux are taken 
from Oort and Rasmusson (1971). Their definitions of stationary and transient 
eddy fluxes have been adhered to In computing the model statistics. 

The apparent maximum in the lower 'troposphere between 30N and 38N is an 
artifact of the reduction of Che fields below the Himalayas. 


- 65 - 


large positive heat flux in the model penetrates higher than is observed, and 
the upper tropospheric maximum is unrealistically weak in the model. 

The stationary eddy heat flux in July is very weak for both the model 
(Figure 5.35b) and the observations (Figure 5.35a). However, the model results 
show several features that are not observed, including maxima of negative flux 
in the subtropical lower troposphere, a maximum of positive flux at high lati- 
tudes in the lower troposphere, and a strong maximum of postive flux at the 
mid-latitude tropopause. The transient eddy heat flux in July, shown for the 
observations (model) in Figure 5.36a (5.36b), is much stronger. The model 
reproduces well the observed broad maximum of positive flux in mid-latitudes, 
although the peak values in the model seem slightly tDO large. Further, the 
maximum of positive flux near the tropopause in mid— latitudes is very weak 
in the model results. 

Local covariance analysis. The observed local transient northward 
heat flux a; 850 mb in winter, as reported by Blackmon et al. (1977) is shown 
in Figure 5.37a. This quantity is given by the covariance of northward velocity 
with temperature, where each field has had the annual cycle removed in the 
manner described in subsection c. The Atlantic and Pacific storm tracks corres- 
pond to regions of high (positive) covariance, and other (positive) maxima 
occur over the Gulf of Alaska, the east coast of Iceland, and southern Russia. 
Minima (areas of relatively lar^je negative flux) are seen over western Canada 
and northwestern Russia. The simulated northern hemispheric winter covariance 
of transient heat flux at 850 mb (Figure 5.37b) shows some similarity to the 
observations. Both major oceanic storm tracks stand out clearly, although the 
heat flux associated with the one in the Atlantic is somewhat too small. The 
maxima near the Gulf of Alaska and in southern Russia are present, although 
their magnitudes, orientations and precise positions do not agree with the 
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observations. The observed minimum over western Canada Is very weak in the 
model results. 

The observed low-pass and band-pass transient 850 mb heat fluxes (obtained 
from the covariances of the filtered v-wind and temperature fields) are shown 
in Figures (5.38a) and (5.39a). The band-pass heat flux is dominated by the 
major oceanic storm tracks, while the low-pass results show the other major 
features that were present in Figure 5.37a. Tvo features that appear as part 
of the storm tracks in Figure 5.37a, but which are seen to correspond to lower 
frequency phene *ena in Figure 5.38a, are the maxima of heat flux over the 
eastern United States and over southeast Asia. These two features are present 
in the simulated 850 mb low-pass heat flux (Figure 5.38b), whose resemblance to 
the tot-il simulated transient heat flux (Figure 5.37b) is close but not as 
close as in the observations. The simulated band-pass heat flux (Figure 5.39b) 
is dominated by the major oceanic storm tracks, which agree very well with the 
observations (Figure 5.39a) both with regard to position and Intensity. 

The observed local transient heat flux in winter at 300 mb and 200 mb 
(Figures 5.40a and 5.41a) present an interesting contrast with one another. 

(Lau, 1978). There are maxima in the northward heat flux over the western 
oceans at both levels, but the heat fluxes over western North America and 
western Europe show an abrupt change in going from large and equatorward at 
300 rob to large and poleward at 200 mb. The simulated winter transient heat 
fluxes at 300 mb and 200 mb (Figures 5.40b and 5.41b) show some evidence of 
this sharp transition. The 200 mb heat flux resembles the observations to the 
extent that it is maximum over the western oceans and western continents, but 
the peak values of the heat flux (except over North America) are unrealistically 
low. In addition, the model results indicate an area of very strong equatorward 
flux in the eastern Pacific and one of more moderate equatorward flux further 
west. Neither of these are observed. The simulated 300 mb heat flux reaches 
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a maximum of negative valuee over western North America with a peak magnitude 
that agrees with the observations. Thus, in this region the model reproduces 
the change of sign of the heat flux. The 300 mb heat flux reproduces the 
observed (positive) maxima off the East coasts of North America and Asia, 
although they are too weak. The observed large, equatorward flux over Western 
Europe at 300 mb is not present in the simulated results. 

Wavenumber-frequency cospectra. The behavior of the transient heat 
flux can be farther described by its wavenumber-frequency cospectrum, which 
expresses the longitude- and time-covariance of v and T in terms of contribu- 
tions from different wavenumbers and frequencies. The model cospectra reported 
here were calculated in a manner entirely analogous to the wavenumber-frequency 
variance spectra described in subsection c. The obse.”ed cospectra are taken 
from Pratt (1975)^^. Figures 5.42a and 5.43a show the observed winter heat 
flux cospectra at 50N at 850 mb and 200 mb, respectively. The corresponding 
cospectra of the GCM are given in Figures 5.42b and 5.43b. Both the observed 
and simulated 850 mb heat flux at 50N show a distinct propagation band, with 
the maximum covariance for wavenumbers 4-5 and for low frequencies. The model 
cospectrum, In general, looks quite realistic, the main exception being the 
lack of covariance near the low frequency end of the propagation band. 

In contrast, the GCM heat flux cospectrum at 200 mb, 50N, looks quite 
unrealistic. The observed cospectrum again shows a well-defined propagation 
band, with a maximum covariance at wavenumber 4, at low frequencies. There is 
also a slight secondary maximum at wavenumbers 7-8 and frequencies near .30, 
space and time scales which are associated with cyclones. The propagation band 
so evident In the observations is not present In the GCM cospectrura, which is 
dominated by very large values at wavenumber 3 near frequencies of 0.10 cpd. 

13 The methodology used by Pratt (1975) is similar to that described In footnote 9. 
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In addition, the low frequency cospectrum for wavenumber 2 is sizeable and 
negative, indicating a relatively large southward heat transport by these 
scales of motion. It Is worth noting that the weak secondary maximum in the 
synoptic-scale region is evident in the model results. 

The summer heat flux cospectrum at 850 mb, SON (Figure 5.44a) also indi- 
cates the presence of a propagation band, but with its low frequency end shifted 
to higher wavenumbers compared to the corresponding winter cospectrum. The 
magnitudes of the heat flux are of course also much less than in winter. This 
decrease in magnitude and shift to higher wavenumber is captured by the GCM 
summer cospectrum at 850 mb, SON (Figure 5.44b), which agrees well qualitatively 
with the observations. There are some differences in detail, with the observa- 
tions showing low frequency peaks at wavenumbers 4 and 6, whereas maxima in the 
simulated cospectrum occur at wavenumbers 6 and 7. 

The propagation band of the observed 200 mb cospectrum at SON (Figure 5.45a) 
is not noticeably weaker than its winter counterpart, although it is shifted 
towards higher wavenumbers. The simulated spectrum (Figure 5.45b) is similar 
to the observed only in that there is some indication of a propagation band 
present. The GCM cospectra show a large drop in peak magnitudes in going from 
winter to summer, and the maximum values in the summer cospectra appear at very 
low frequencies in wavenumbers 2 and 5. These features do not correspond to 
the observations. 

e. Northward eddy transport of momentum-stationary and transient eddy 
momentum flux. The time mean momentum flux due to the stationary eddies in 
February is shown (as a function of latitude and pressure level) for the obser- 
vations^^ in Figure 5.46a, and for the GCM in Figure 5.46b. The observations 

The observed results for stationary and transient momentum flux are taken 
from Oort and Rasmusson (1971). Their definitions of stationary and transient 
eddy fluxes have been adhered to in computing the model statistics. 
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are dominated by large positive (northward) values of the flux in the upper 
troposphere at raid-latitudes, and large negtive (southward) values farther 
north. The peak northward flux in the model upper troposphere in raid-latitudes 
is larger than observed, while the observed region of large southward flux is 
almost completely missing in the GCM results. 

The observed time mean momentum flux due to transient eddies for February 
is shown In Figure 5.47a. It is dominated by a large peak in northward flux 
slightly north of SON at the 200 mb level. The only other feature of note is 
the relatively small peak of negative (southward) flux at 200 mb at the equator. 
The model results (Figure 5.47b) also show a large maximum of positive flux, 
and its position is correct. However, it is too weak by a factor of two. A 
small maximum of southward flux at 200 mb over the equator is also present in 
the model, but it is too weak. A feature present in the GCM results which is 
not observed Is the region of negative flux in the upper troposphere between 
54 and 70N. 

The behavior of the GCM’s time mean momentum flux due to the stationary 
eddies for July is quite unrealistic (Figures 5.48a-b). The observations show 
a peak of northward flux at 200 mb at about 30N, while in the GCM results this 
peak is completely missing. The (small) maxima of negative flux in the model 
upper troposphere do not have their counterpart in the observations. In con- 
trast to the stationary eddy flux of momentum, the model transient eddy flux of 
momentum behaves more realistically. The observations (Figure 5.49a) show 
three important features: a large maximum of positive flux at 200 mb at 43N, a 

smaller (positive) maximum at 200 mb just south of the equator, and a broad 
upper tropospheric maximum of negative flux north of 60N. The GCM results 
(Figure 5.49b) show the same basic features, located approximately in the cor- 
rect positions. However, the positive maxima are too small, and the negative 
maxiraun too large in magnitude. 
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Local covariance analysis. The observed local transient northward 


flux of momentum at 250 mb In winter, as reported by Lau et al» (1978), is 
shotm In Figure 5.50a. This quantity Is given by the covariance of northward 
velocity with eastward velocity. The covariance of the band**pass filtered 
fields is given In Figure 5.51a. The main features of these two observed 
fields are similar, although the magnitude of the band~pass covariance is con- 
siderably less than that of the total covariance. One feature that Is prominent 
in both maps is the dipole structure that appears over the western and central 
Atlantic with poleward fluxes south of about 45N and equatorward fluxes north 
of this. This feature Implies a large convergence of momentum flux into the 
storm track region (defined by the maximum RMS of the 500 mb band-pass height 
field), which Is Indicated In the figures by the dashed arrows. A similar 
feature appears in the Pacific. Another Important feature is the large north- 
ward flux apparent over the western parts of continents, particularly the 
United States. The GCM local winter covariance at 200 mb is shown in Figure 
5.50b, and the band-pass local covariance in Figure 5.51b. The simulated map 
of local covariance shows dipole patterns over North America and the Eastern 
Pacific, off the east coast of Asia, and over Western Europe. Clearly the sim- 
ulated convergence of momentum flux is not realistic. The Pacific and Atlantic 
storm tracks seem to be marked by almost no convergence of momentum flux, in 
contrast to the observations. The GCM band-pass covariance has dipole-like 
patterns in the Atlantic (although it is south of the storm track) and in the 
eastern Pacific, but they are too weak. The maximum over western North America 
is present in the GCM results but it is too weak. 

The observed winter covariance of westward and northward velocities at 
500 mb (taken from Blackmon et al. , 1977) is given in Figure 5.52a, and the band- 
pass variance in Figure 5.53a. The map of total covariance shows a general pat- 
tern of momentum flux convergence into mid-latitudes, although the dipole-like 
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structures evident st 250 mb are not clear here. Features tfhlch the model 
covariances (Figure 5. 52b) capture are the peaks of poleward flux over the 
extreme western oceans and western continents (although the magnitude of the 
model flux Is generally too small), and the tongue of equatorward flux near 
Greenland. 

The observed winter band-pass variance at 500 mb Is similar to the 250 mb 
covariance, with distinct dlpole-llke patterns over the western oceans, and a 
large poleward flux over Che western United States. The GCM band-pass covari- 
ance (Figure 5.53b) at 500 mb is much more realistic than at 200 mb. The dipole 
pattern observed In Che Atlantic Is evident In the model results, although It 
Is too extensive. In the Pacific, the observed dipole pattern Is also evident 
in Che model results, but It appears to be too diffuse. 

Wavenumber-frequency cospectra. The transient momentum flux can be 
further examined by means of wavenumber-frequency cospecCra, exactly analagous 
to Che heat flux cospectra described in subsection d. The observed winter 
wavenumber-frequency cospectrum of northward and eastward velocities at 300 mb, 
35N (taken from Pratt, 1975) Is shown in Figure 5.54a, and the GCM cospectrum 
(at 34N) is shown in Figure 5.54b, Both the observed and the simulated cospectra 
show a clear propagation band. The observed cospectrum peaks at wavenumber 5, 
at a frequency of 0,10 cpd, whereas the model results show two low frequency 
peaks, at wavenumbers 2 and 5. The presence of a double peak makes the slope 
of the model propagation band ambiguous at the low frequency, low wavenumber 
end. The observed region of negative (equatorward) flux at low frequencies and 
wavenumbers 8 and 9 Is present in the model, although the choice of contours in 
Figure 5.54b obscures this somewhat. The secondary peak at wavenumber 7 and 
frequencies near 0.175 cpd. does not appear in the observations. 
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Balance of transient »ddy flux converg»nce and ageoitrophlc accelera- 
tion In the local noaentua budget . In the zonally-averaged budget of zonal 
momentum, the contribution of the Corlolla acceleration due to the time-mean 
meridional ageoatrophlc wind la very weak. However, thla Is no longer true in 
the framework of a local budget (Lau, 1978). In fact, locally the advectlon of 
zonal momentum by the time-mean circulation is basically balanced by the ageo- 
strophlc Coriolis acceleration, the horizontal convergence of transient momentum 
flux playing a less important role. This (two-dimensional) convergence Is cal- 
culated from observations for winter at 250 mb in Lau (1978), and Is presented 
here in Figure 5.55a. This should be contrasted to the Coriolis acceleration 
accompanying the time-mean meridional ageostrophlc flow shown in Figure 5.56a. 

In the vicinity of the two major jet streams (indicated in Figure 5.56a by the 
solid arrows) the ageostrophlc acceleration dominates the transient flux contri- 
bution, causing a local acceleration in the jet entry regions and a local deccel- 
eration in the jet exit regions. In the area of the North Amerlcan-Atlantlc jet 
stream, the transient fluxes work to oppose these accelerations. Their effect 
in the region of the Japan jet is not as systematic. It should be pointed out 
that the dominance of the ageostrophlc term is seen also over large portions 
of Europe and Asia. 

The corresponding GCM results for the convergence of transient momentum 
flux and for the ageostrophlc Coriolis acceleration at 285 mb for the month of 
February are shown in Figures 5.55b and 5.56b. Figure 5.56b also shows the 
model's jet streams. It is apparent that the Pacific jet is correctly located, 
but that the Atlantic jet is not oriented towards the northeast as in the 
observations. In comparing the model results to those that are observed, it is 
the relationship of the terms in the momentum budget to the jet streams that 
will be considered. (Note also that the observed results are obtained from 11 
winter seasons, each winter being defined as the 120 day period starting on 
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15 November, whereas the model results represent an average over one month 
only.) The ageostrophic Coriolis acceleration in the GCM acts very realisti- 
cally in relation to the Jets, causing an acceleration over eastern Asia and 
off the southern coast of the United States, and causing a deccele ration in the 
jet exit regions. The magnitude of these accelerations is quite realistic for 
the Atlantic jet, but is too weak in the Pacific. The horizontal convergence 
of transient momentum flux is realistic in the sense that it is of the observed 
magnitude, approximately a factor of two less than the ageostrophic accelera- 
tions. Otherwise the patterns show differences. The observed convergence 
(acceleration) in the Atlantic is replaced in Che model with a divergence, while 
the observed convergence over western North America appears in the eastern 
Pacific in the GCM. The observed divergence over Japan is realistically simu- 
lated, but the observed divergence downstream is replaced in the model by a 
region of convergence. Overall, the model's transient momentum flux convergence 
does not have a consistent relationship to the positions of the jets. 
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FlKure 5.1 Energy cycles for the observed and simulated winter hemispheres, 

computed in the mixed space-time domain. is the zonal-mean 

available potential energy, Pg the eddy availa^'le potential 
energy. Kg the eddy kinetic energy, and Kftf the zonal-mean kinetic 
energy. The energies are quoted in units o^ 10^ joules/m^, 
the transformations and sources and sinks in units of watts/m^. 


j (a) Observed Northern Hemisphere, January, (b) GCM Northern 
n A Hemisphere, February, (c) GCM Southern Hemisphere, July. (Values 

PACE li^TK)lMU.Y''kSfAr‘'’"'' .ppro«i..te.) 
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Figure 5.2 Energy cycles for the observed and simulated summer hemispheres, 

computed as In Figure 5.1, (a) Observed Northern Hemisphere, July, 
(b) GCM Northern Hemisphere, July, (c) GCM Southern Hemisphere, 
Ftbiuary. (Values In parentheses are very approximate) 
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Figure 5.3b GCM niorldlonal mean energy transports for February, units of 
10^^ calories/day , given as in Figure 5.3a. V/here transports 
are not shown, they are extremely small. 
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Figure 5.3c Observed meridional eddy energy transports for February, units 

of 10^^ calorl es/day, given as in Figure 5.3a. Where transports 
are not shown, they are extremely small. 
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Figure 5.3d GCM meridional eddy energy transports fo* lr*oruary, units of 

calorles/day , given as in Figure 5.3a. Whsre transports 
are not shown, they are extremely small. 
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Figure 5.4c Observed meridional eddy energy transports for July, units of 
in^^ calories/day , given as in Figure 5.3a. Wliere transports 
are not shown, they are extremely small. 
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Figure 5.4d GCM meridional eddy energy transports for July, units of lO^^ 
calories/day , given as in Figure 5.3a. Where transports are 
not shown, they are extremely small. 
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Figure 5.5bl GCM winter stationary planetary wave variance of the geopotential 
height in the Northern Hemisphere. Units of 10^ m2. 
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Figure 5.6bl r.Ot winter amplitude of the geopotential height, stationar; wave 
number 1, in the Northern Hemisphere, units of 10^ m. 
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Figure 5.6dl GCM winter amplitude of the geopotential height, stationary 
wavenumber 2, in the Northern Hemisphere, units of 10^ m. 
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Figure 5.<bel Observed winter amplitude of the geopotentlal height, stationary 
wavenumber 3, In the Northern Hemisphere, units of 10 *• 
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Figure 5.6e2 Observed winter phase of the geopotential height, stationary 
wavenuaber 3, in the Northern Henisphere, units of degrees, 
phase is relative, so that the longitude of the ridge is the 
phase divided by the wavenuaber. 
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Figure 5.7a Observed summer stationary planetary wave variance of the 

geopotential height in the Northern Hemisphere. Units of 10-^ 
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Figure 5.7bl GCM summer stationary planetary wave variance of the geopotentlal 
height In the Northern Hemisphere. Units of 10^ m^. 
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Figure 5.7b2 GCM sumaer stationary planetary wave variance of the geopotential 
height in the Southern Hemisphere. Units of 10^ 






Figure 3«Sa2 Observed suonner phase of the geopotential height, stationary 
wavenumber I, In the Northern Hemisphere, units of degrees, 
phase Is relative, so that the longitude of the ridge Is the 
phase divided by the wavenumbet. 





Figure 5.8bl GCM summer amplitude of the geopotential height, stationary 
wavenumber 1, in the Northern Hemisphere, units of 10 ra. 
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Figure 5.8b2 GCM Bu-er phaBe of the geopotantlal height, at«lonary -a/enuaber 
I In the Northern Hemisphere, units of 10 m. The phase is 
relative, so that the longitude of the ridge Is the phase divided 

by the wavenumber. 
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Figure 5.8cl Observed Sumer amplitude of the geopotential height, stationary 
wavenumber 2, in the Northern Hemisphere, units of 10 m. 
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Figure 5.8c2 Observed summer phase of the geopotential height, stationary 
wavenumber 2, in the Northern Hemisphere, units of degrees, 
phase is relative, so that the longitude of the ridge is the 
phase divided by the wavenumber. 




Figure 5.8dl GCM summer amplitude of the geopotential height, stationary 
wavenumber 2, in the Northern Hemisphere, units of 10 m. 
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Figure 5.8d2 GCM summer phase of the geoptentlal height, stationary wavenumber 
2, In the Northern Hemisphere, units of degrees. The phase Is 
relative, so that the longitude of the ridge Is the phase divided 
by the wavenumber. 
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Figure 5.8el Observed summer amplitude of the geopotential height, stationary 
wavenumber 3. in the Northern Hemisphere, units of 10 m. 
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Figure 5.8e2 Observed summer phase of the geopotential height^ stationary 
wavenumber 3, in the Northern Hemisphere, units of degrees, 
phase is relative, so that the longitude of the ridge is the 
phase divided by the wavenumber. 
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Figure 5.8f2 GCM summer phase of the geopotential height, stationary wavenumber 
3, in the Northern Hemisphere, units of degrees. The phase is 
relative, so that the longitude of the ridge is the phase divided 
by the wavenumber. 
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Figure 5.9b GCM RMS deviation of the winter 500 mb geopotential height field 
(Northern Hemisphere), units of m. Contour interval is 10 ■. 
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Figure 3.10a Observed RMS deviation of the wintet 300 mb lotr-paas geopotential 
height field (Northern Hemisphere), units of m. Contour interval 



Figure 3.10b GCM RMS deviation of Che winter 500 ab low-pass geopotential height 
field (Northern Henlsphere) , units of ■. Contour Interval Is 10 a. 
(Several contours have jeen omitted In the Immediate vicinity of 
very sharp maxima or minima.) 
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Figure 5.11a Observed RMS deviation of the winter 500 mb band-pass geopotential 
height field (Northern Hemisphere), units of m. Contour Interval 
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Figure 5.11b GCM RMS deviation of the winter 500 mb band-pass geopotential height 
field (Northern Hemisphere), units of m. Contour interval is 5 m. 
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Figure 5.12a Observed RMS deviation of the sumaKr 500 mb geopotential height 
field (Northern Hemisphere), units of m. Contour Interval is 10 
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Figure 3«13a Observed RMS deviation of the sununer band~pas8 500 mb geopotential 
height field (Northern Hemisphere), units of m. Contour interval 
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Figure 5.13b GCM RMS deviation of the summer band-pass 500 mb geopotential height 
field (Northern Hemisphere), units of m. Contour Interval is 5 m. 
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Figure 5.14b GCM RMS deviation of the winter 850 mb temperature field (Northern 
Hemisphere), units of degrees K. Contour interval is 1 degree. 
(Several contours have been omitted in the immediate vicinity of 
very sharp maxima or minima.) 
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Figure 5.15b GCM RMS deviation of the winter low-pass 850 mb temperature field 
(Northern Hemisphere), units of degrees. Contour interval Is 
0.5 degrees. (Several contours have been omitted in the imnediate 
vicinity of very sharp maxima or minima.) 
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Figure 3.16b GCM RMS deviation of the winter band-pass 850 mb temperature field 
(Northern Hemisphere), units of degrees K. Contour interval is 
0.3 degrees. (Several contours have been omitted in the imaediate 
vicinity of very sharp maxima or minima.) 
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Figure 5.17b GCM wavenumber-frequency spectral density of winter 500 mb 

geopotential height at 50N, in units of 10^ m2/cpd. Otherwise 
as in Figure 5.17a. 
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Figure 5.18a Observed wavenumber-frequency spectral density of winter 500 rat 
geopotential height at 30N, in units of 10 ra /cpd. Otherwise 
in Figure 5.17a. 
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ZONAL WAVENUMBER 



Figure 5.19b GCM wavenumber-frequency spectral density of summer 500 mb 

geopotential height at SON, in units of lO^ m2/cpd. Otherwise 
as in Figure 5.17a. 
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Figure 5.20a Observed wavenumber-frequency spectral density of summer 500 mb 
geopotentlal height at SON, in units of 10 in^ cpd* Otherwise 
as in Figure 5.17a. 
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Figure 5.20b GCM wavenumber-frequency spectral density of summer 500 mb 

geopotential height at 30N, in units of 10^ m^/cpd. Otherwise 
as in Figure 5.17a. 
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Figure 5.21a Observed wavenumber-frequency net propagation tendency of winter 
500 mb geopotential height at SON, dimensionless. 
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Figure 5.21b GCM wavenunber-f requency net propagation tendency of winter 500 
mb geopotential height at SON, dimensionless. 



(Qil3) A3N3003UJ 


148 


Figure 5.22a Observed wavenumber-frequency net propagation tendency of winter 
500 mb temperature at 50N, dimensionless. 
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Figure 5.23a Observed wavenumber-frequency spectral density of winter 300 mb 
u-wind at 35N, in units of m2s“2/cpd. Otherwise as in Figure 
5.17a. 
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Figure 5.23b GCM wavenumber-frequency spectral density of winter 300 mb u-wind 
at 34N, in units of m^s'^/cpd. Otherwise as in Figure 5.17a. 
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Ficure 5.24a Observed wavenumber-frequency spectral density of winter 300 
v-wind at 35N, in units of m^s"2/cpd. Otherwise as in Figure 
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Figure 5.24b GCM wavenumber-frequency spectral density of winter 300 mb v-wir.H 
at 34N, In units of m2s“2/cpd. Otherwise as in Figure 5.17a. 



15 YEAR AVERAGE : LOW FREQUENCY PLANETARY WAVES 



Figure 5.25a Observed winter geopotential height variance in low-frequency 
planetary waves, in units of loW (Northern Hemisphere). 
Contour interval is .1 x lO^m^, 
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Figure 5.25b GCM winter geopotentlal height variance in low-frequency plane- 
tary waves, in units of lO^ni^ (Northern Hemisphere). Contour 
Interval is . 1 x 10^ m^. 
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Figure 5.26a Observed winter geopotential height variance in low-frequency 
synoptic-scale waves, In units of 10^^ (Northern Hemisphere). 
Contour interval is .4 x lO^ra^. 
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Figure 5.26b GCM winter geopotential height variance in low-frequency synoptic- 
scale waves, in units of lO^a^ (Northern Hemisphere). Contour 
Interval is .4 x 10^ m^. 
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Fiffure 5.30b GCM summer geopotential height variance in low frequency synoptic- 
scale waves, in units of 10^ (Northern Hemisphere). Contour 
interval is .A x 10^ m^. 



200 



166 


Figure 5.31a Observed summer geopotential height variance in medium frequency 
planetary waves, in units of 10^ (Northern Hemisphere). 
Contour Interval is .1 x 10^ m^. 
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Figure 5.31b GCM suraner geopotentlal ’..eight variance in aedlua frequency 
planetary waves, in units of 10 n (Northern Heaisphere)* 
Contour interval is 0.1 x 10^ m^. 
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Figure 5.32a Observed summer geopotential height variance in medium frequency 
synoptic-scale waves, in units of 10^ m^ (Northern Hemisphere), 
Contour interval is .1 x 10^ 
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Figure 5.32b GCM summer geopotential height variance In medium frequency 

synoptic-scale waves, In units of 10^ m^ (Northern Hemisphere) 
Contour Interval Is .2 x 10^ ra^. 
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Figure 5.33b GCM February stationary eddy heat flux, In units of “K m/sec 
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Observed February transient eddy heat flux. In units of *K m/sec. 



GCM February transient eddy heat 



flux, in units 



o 


(D 

OD 



OOOOOOOOO 

OOOOoOOOO 


<u 

T3 

3 



(quj) djnssdJcj 
% 


175 


Figure 5.35b GCM July stationary heat flux. In units of “K m/sec. 
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Figure S.36b GCM July transient eddy heat flux, in units of *K a/sec. 




Figure 5.37a Observed winter local transient heat flux at 850 mb, in units of 
"K m/sec. Contour Interval is 10“K m/sec (Northern Hemisphere). 



Figure 5.37b GCM winter local transient heat flux at 850 mb, in units of 

"K a/sec. Contour interval is 10“K m/sec (Northern Hemisphere). 
(Several contours have been omitted in the immediate vicinity of 
very sharp maxima or minima.) 
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Figure 5.39b GCM winter local band-pass transient heat flux at 850 mb, in 

units of “K m/sec. Contour interval is 2®K m/sec (Northern Hemi- 
sphere). 
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Figure 5.40b GCM winter local transient heat flux at 300 mb, in units of 

“K m/sec. Contour interval Is 5“K m/sec (Northern Hemisphere). 
(Several contours have been omitted in the immediate vicinity oJ 
very sharp maxima or minima.) 
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Figure 5.42b GCM wavenumber-frequency cospectral density of winter 850 mb 
heat flux at 50N, in units of “K m/sec/cpd. Otherwise as in 
Figure 5.42a. 
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Figure 5.43b GCM wavenumber-frequency cospectral density of winter :^00 mb 
heat flux at 50N, In units of »K m/sec/cpd. Otherwise as in 
Flgur® 5.42a, 
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Figure 5.44b GCM wavenumber-frequency cospectral density of summer 850 mb 
heat flux at 50N, In units of °K m/sec/cpd. Otherwise as In 
Figure 5.42a. 
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Figure 5.45b GCM wavenuraber-f lequency cospectral ^-^nsity of sumner 200 mb 
heat flux at 50N, In units of “K m/sec/cpd. Otherwise as in 
Figure 5.42a. 
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Figure 5.46b GCM February stationary eddy momentum flux, In 
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Figure 5.A7b GCM February transient eddy momentum flux, in units of (m/sec)2 
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Observed July stationary eddy Domentum flux, in units of (m/sec)2. 
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Observed July transient eddy momentum flux, in units of (m/sec)2. 
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transient eddy aonentun flux, in units of (m/sec) 






205 

I 






207 


Figure 5.51b GCM winter local band-pass transient aoiaentua flux at 200 ab In 
units of (a/sec)2. Contour Interval Is 5 (a/sec)^ (Northern 
Healsphere) . 





Figure 5.52b GCM winter local transient momentum flux at 500 mb in units of 

(m/sec)2. Contour interval is 20 (m/sec)^ (Northern Hemisphere) 
(Several contours have been omitted in the Immediate vicinity of 
very sharp maxima or minima.) 
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Figure 5.53b GCM winter local band-pass transient momentum flux at 500 mb in 
units of (m/sec)2. Contour Interval is 4 (m/sec)^ (Northern 
Hemisphere). (Several contours have been omitted in the immediate 
vicinity of very sharp maxima or minima.) 
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Figure 5.54a Observed wavenumber-frequency cospectral density of winter 300 mb 
momentum flux at 35N, in units of (m/sec)2/cpd. Otherwise as 
in Figure 5.42a. 
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Figure 5.54b GCM wavenumber-frequency cospectral density of winter 300 mb 
momentum, flux at 34N, In units of (m/sec)Vcpd. Otherwise as 
In Figure 5.42a. 
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Figure 5.35a The observed local horizontal convergence of transient flux of 
westerly noaentuffl for Northern Healsphere winter at 250 mb. 
Plotted In units of 10“^ m/sec^, with a contour Interval of 50 
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Figure 5.55b The GCM local horizontal convergence of transient flux of westerly 
momentum for the Northern Hemisphere during February at 285 mb. 
Plotted in units of 10“^ m/sec^, with a contour interval of 100. 
(Several contours have been omitted in the immediate vicinity of 
very sharp maxima or minima.) 
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6. Simulation of dlabatlc and frictional processes 


a. Summary of the global sensible heat and moisture budgets. In both sim- 
ulations all fluxes of heat and moisture across the top and bottom boundaries 
of the atmosphere were accumulated at every grid point. Global sums were car- 
ried out to obtain the average global values of each boundary flux. Figs. 6.1a-b 
and Tables 6.1 and 6.2 summarize the global sensible heat and moisture budgets 
for the February and July simulations, giving the globally averaged surface 
sensible heat flux, evaporation, their ratio (termed the global Rowen ratio), 

the precipitation rate, and the planetary albedo. The observed surface fluxes 
are these of Budyko (1963) as given by Schutz and Gates (1971, 1972), the 
observed precipitation rates are from Jaeger (1976), and the observed planetary 
albedoes are from Raschke, et al. (1973). The simulated global mean evaporation 
and precipitation rates agree well with observations, and are Improved relative 
to earlier versions of the model. The simulated surface sensible heat flux Is 
also significantly Improved, but It continues to be excessive. Finally, the 
simulated planetary albedoes are too high In both February and July, and are 
Inferior to these obtained with earlier versions of the model. 

The following subsections explore these and other results In more detail. 

b. Surface fluxes . The simulated and observed (Schutz and Gates, 1971, 
1972) January/February and July geographical distributions and zonal means of 
the evaporation are shown In Figs. 6.2a-f. In the February simulation, the 
evaporation maxima off the east coasts of Asia and North America have their 
observed magnitudes, although their shapes are not well-simulated and the 
secondary maximum south of Greenland Is missed. The observed equatorial minimum 
Is simulated In the Eastern Pacific, but Is missed in Melanesia and the Indian 
Ocean. The zonal mean does not show the observed sharp equatorial minimum. 
Observed maxima In the tropical South Atlantic, Indian, and South Eastern 


-'*’5. / r' 

PACE mUHTIONAUV BUANA 


Pacific oceans are either missing or weak and malformed. 

Similarly, the simulated and observed (Schutz and Gates, 1971, 1972) heat 
flux distributions and zonal means for January/February and July are shown In 
Figs. 6.3a-f. The February simulated surface sensible heat flux has maxima of 
the observed Intensities off the east coasts of Asia and North America, but the 
Atlantic maximum falls to extend to the west coast of Scandinavia; it hugs the 
coast of North America more closely than observed. 

In the July simulation, the evaporation maxima off the coasts of Brazil 
and Madagascar are correctly simulated, as are the minima over Australia and 
North Africa. To some extent, the minima over such desert areas must be cor- 
rectly simulated, since the ground moisture distribution was Initialized from 
observations and has a memory of several weeks. The general pattern of oceanic 
evaporation Is not very satisfactory. The simulated maximum In the tropical 
Western Pacific Is too strong, and does not show the observed relative minimum 
at the equator. The equatorial minimum Is also missed In the central Pacific. 

The zonal mean of the simulated surface sensible heat flux Is much stronger 
then observed. The simulated relative minimum at 10* north Is about 20 W m"^. 

The largest zonal mean values occur at 30* north. In both the simulation and 
the observations, but the simulated maximum Is too strong. 

Figs. 6.4b and d show the simulated global surface stress distribution for 
February and July, and Figs. 6.4a and c show the corresponding observed stresses, 
over the oceans only, as derived by Mlntz (1979) from the data of Hellerman 
(1967). The simulation Is reasonably successful In reproducing the observed 
stresses, except In the "roaring fifties" around Antarctica, where the simulated 
stresses are weak and disorganized, reflecting the model's failure to simulate 
the dynamics of that region. 

c. Precipitation . The simulated and observed (Jaeger, 1976) February and 
July distributions of precipitation are shown In Fig. 6.5a-d. The simulated 


precipitation field la In excellent agreement with observations. Almost all 
Important features are well-simulated In both position and Intensity. Perhaps 
the single major discrepancy is the somewhat excessive (2 mm/day too large) 
February precipitation over the Tibetan plateau. In the July simulation the 
major monsoons are captured. The observed maximum over eastern North America 
is somewhat overpredicted. 

Figs. 6.5e and f show the partitioning of the simulated total precipitation 
between large-scale-saturation and convection, In the zonal mean, for February 
and July. In both seasons, cu.^ulus precipitation dominates In the tropics, 
while large-scale precipitation dominates In higher latitudes. 

d. Cloudiness and albedo. The simulated and observed February and July 
distributions of total cloudiness are showm In Fig. 6.6a-d. It Is apparent chat 
the simulated cloudiness Is excessive. This Is confirmed by the fact that the 
simulated planetary albedoes are too high (see again Tables 6.1-2): the simu- 

lated February and July planetary albedoes are 0.390 and 0.380, while the corre- 
sponding observed values (Raschke et al. , 1973) are 0.283 and 0.281, respectively. 
Much of the excessive cloudiness is due to widespread "supersaturatlon” in the 
lowest model layer. Moeng and Randall (1982) show that even If cumulus cloud- 
iness Is completely neglected, the simulated albedo Is practically unchanged. 
Apparently the boundary layer and cumulus parameterlzatlons of the model do not 
allow realistic transport of moisture up out of the lowest model layer (Randall, 
1982). 

The model correctly simulates Che seasonal Increase In zonal ly-averaged 
cloudiness from summer Co winter In the high latitudes of each hemisphere, 
except that the observed Increase In Arctic cloudiness from February to July 
Is missed. The cloudiness of the ITCZ Is reasonably well-simulated, particul- 
arly for February, and the observed subtropical minima are also captured. 
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The simulated longitudinal distribution of cloudiness Is much less satis- 
factory. The observed subtropical July maxima off the west coasts of North 
America, South America, and South Africa are not only missed; they ere replaced 
by minima . The observed maxima are known to be associated with stratus and 
stratocumulus clouds within the planetary boundary layer. The fact that the 
model has not simulated these maxima Indicates that an Improved planetary 
boundary layer parameterization Is needed (Randall, 1982). 

e. Total dlabatlc heating. Fig. 6.7a shows the "observed" zonally aver- 
aged total dlabatlc heating reported by Newell et al. (1972), and Figs. 6.7b 
and c show the corresponding simulations for February and July. The simulations 
produce deep equatorial maxima associated with the convection of the ITCZ, and 
prominent low-level winter-hemisphere maxima due to boundary layer heating. 

The middle and upper troposphere are radlatlvely cooled. In each season, cool- 
ing predominates above 800 mb In the subtropics of the winter hemisphere. This 
radiative cooling Is balanced by the adiabatic warming In the descending branch 
of the Hadley circulation (see Section 7). The low-level equatorial heating 
reported by Newell et al. is negative, but It Is positive In the model. 
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Table 6.2a February energy balance (W m”^). 



Analysis of heat budget fields 

Items 

S. hemisphere 

N. hemisphere 

Global 


Mean 

Mean 

Mean 

Top Boundary of Atmosphere 




Incoming short wave flux 

440.18 

273.14 

356.83 

Outgoing short wave flux 

178.59 

100.09 

139.34 

Outgoing long wave flux 

196.12 

194.24 

195.18 

Within the Atmosphere 




S.U. absorbed 

65.18 

40.80 

52.99 

L.K. release by precip. 

99.67 


89.19 

S.H. increase 

1.12 

2.16 

1.64 

KE Increase 

0.05 

-0.06 

-0.005 

LE increase 

.70 

.66 

.68 

Bottom Boundary of Atmosphere 




S.W. absorbed 

196.40 

132.60 

164.50 

L.W. outgoing 

27.78 

35.49 

31.63 

S.H. flux outgoing 

23.53 

31.55 

27.54 

L.H. flux outgoing 

88.78 

87.54 

88.16 

Net Energy Absorbed 




Land 

— — 

- _ 

-1.30 

Ocean 

- 

- 

25.08 

Lower boundary 

56.32 

-21.98 

17.17 

Top boundary 

65.47 

-20.85 

22.31 

Planetary Albedo 

0.41 

0.37 

0.39 
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Table 6.2b July enerj^y balance (W m"2) 



Analysis of heat budget fields 

Items 

S. hemisphere 

N. hemisphere 

Global 


Mean 

Mean 

Mean 

Top Boundary of Atmosphere 

Incoming short wave flux 
Outgoing short wave flux 
Outgoing long wave flux 

214.99 

76.14 

199.78 

458. 18 
180.80 
210.02 

336.58 

128.47 

204.90 

Within the Atmosphere 




S.U. absorbed 

33.30 

69.30 

51.30 

L.H. release by preclp. 

84.30 

103.76 

94.03 

S.H. Increase 

-1.24 

6.15 

2.45 

KE Increase 

0.098 

0.056 

0.021 

LE Increase 

0.044 

2.729 

___ _| 

1.387 

Bottom Boundary of Atmosphere 




S.W. absorbed 

105.54 

208.09 

156.81 

L.W. outgoing 

32.76 

35.75 

34.26 

S.M. flux outgoing 

31.54 

33.42 

32.48 

L.H. flux outgoing 

95.72 

91.37 

93.54 

Net Energy Absorbed 




Land 

- 

- 

3.60 

Ocean 

- 

- 

-6.48 

Lower boundary 

-54.48 

47.54 

-3.47 

Top boundary 

60.94 

-67.38 

-3.22 

Planetary Albedo 

0.35 

0.39 

0.38 
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Figure 6.2b February simulated evaporation, in mm dy 
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Figure 6.2e July simulated evaporation, mm dy 
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'Served surface sensible heat flux. In W m taken 
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Figure 6.3c January observed zonal nean surface sensible heat flux taken froa 
Schutz and Gates (1971) (solid line), and February siaulated 
zonal mean surface sensible neat flux, Wb" 2 (dashed line). 
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Figure 6.3f July observed zonal mean surface sensible heat flux taken from 
Schutz and Gates (1972) (solid line), and July simulated zonal 
mpan surface sensible heat flux, W m“2. 




Figure 6.4a February observed surface stress, dynes cm~2^ taken from Mintz 
( 1979) as derived from 'lellerman (1967). 
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Figure 6.4b February simulated surface stress. The arrows indicate direction 
and the contour lines indicate magnitude. The contour interval 
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Figure 6.4c July observed surrace scress, ayn 
as derived from Hellerinan ( 1967). 
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Figure 6.'Sb February simulated precipitation. The contour Interval is 
1 mm dy“^ except in congested regions. 
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Figure 6.5c July observed precipitation, taken from Jaeger (1976). The coa 
tour Interval is 1 mm dy“^. 
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Figure 6.6b February slaulated total cloudiness, in tenths. Regions with 
less than 7/10 cloudiness are shaded. 
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Figure 6.6d July simulated total cloudiness. Regions with less than 7/10 
cloudiness are shaded. 
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Figure 6.7b February simulated zonally averaged total diabatic heating, K day” 
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7. Large-scale overturnings 


a. Mean meridional circulations. The mass transport streamfunctlons of 
the simulated February and July mean meridional circulations have been deter- 
mined by Interpolating the six-hourly simulated meridional winds to Isobarlc 
surfaces, time- and zonally-averaglng the result, and subtracting away Its ver- 
tical mean. The streamfunctlon for the time- and zonally-averaged meridional 
circulation Is then given by 

PS 


fH ((^,p) - - (1/g) ! (vl 2^a cos dp' , 

P 


(7.1) 


where v Is the meridional wind, the overbar denotes a time average, [ 1 denotes 
a zonal average, a is the radius of the earth, ^ Is latitude, p is pressure, 
and subscript S denotes the earth's surface. The simulated latitude-height dis- 
tributions of Vjj for February and July are given in Fig. 7. lb and c; the 
observed distributions obtained by Newell et al . (1972) are shown for comparison 
In Fig. 7.1a. In both seasons, but particularly in July, the simulated Hadley 
cells are somewhat weaker than observed. In view of the present model's lack 
of cumulus friction, and the widespread belief (e.g., Helfand, 1979) that 
cumulus friction acts to strengthen the Hadley circulation, this result is to 
be expected if not welcomed. The simulated mid-latitude Ferrel cells are also 
weaker than observed. 


b. The Walker circulations. Bjerknes (1969) suggested the possible Impor- 
tance of large-scale equatorial east-west overturnings, which he called Walker 
circulations. It Is believed (Cornejo-Garr ido and Stone, 1977) that the Walker 
circulation Is driven by a longitudinal variation of the latent heating. Newell 
et al . ( 1972) suggested that the term "Walker circulation" be generalized to 
Include any longitudinal overturning, tropical or otherwise. 
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We have computed the simulated Walker mass transport streamfunctlons by 
the following procedure: The zonal wind field Is Interpolated to Isobaric 

surfaces, and subsequently time-averaged. At each level and latitude the zonal 
mean of the result Is then subtracted away. Finally, the vertical mean of the 
result Is subtracted away at each latitude and longitude to yield closed circul- 
ations. For convenience, let the residual zonal wind obtained by this procedure 

S 

be denoted by u^ , where the overbar again denotes a time mean, the * denotes a 
departure from the zonal mean, and the tilde denotes a departure from the verti- 
cal mean. The mass transport stream function Is then given by 




(1/g) / a dp' 


1,2 


(7.2) 


where ( denotes an average over latitude from to ^2* 

1 ,2 

For both February and July, we have computed \ for the latitude strips 
30®S-Equator and Equator-30"N; the results are shown in Fig. 7.2b-e. For com- 
parison, Fig. 7.2a shows with the observations reported by Newell et al . (1972), 
for 5®S and 5®N. The most striking finding In the simulations Is the very 
large, well-organized Northern Hemispheric circulation In February. The rising 
branch of this Walker cell Is near 140® west longitude, over the eastern North 
Pacific. As shown in Fig. 6.5, there Is an intense precipitation maximum at 
about 10®N at this longitude. The sinking branch of the Walker cell is near 
90® east longitude. In the vicinity of the Siberian high. In July, the simulated 
Northern Hemisphere Walker circulation is much weaker and less well-organized. 

In the Southern Hemisphere winter, there is a modest Walker cell whose 
rising branch is near the spectacular precipitation maximum on the date line, 
and whose sinking branch is near South Africa. The February simulation exhibits 
a weaker, less well-organized Southern Hemispheric Walker circulation. 
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Suiamarlzlngt the model has produced vigorous winter-hemisphere Walker cir- 
culations in both February and July. The cells have their rising branches over 
oceanic regions of strong latent heat release, and their sinking branches over 
land. The summer-hemisphere circulations are much weaker in both February and 
July. 
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Figure 7.2c July siraalated stream function for the Walker circulations 
the Equator to J0*N, In units of 10^ kg s~*. 










8. Summary and Conclusions 


The GLAS climate model simulation of the mean dynamical fields Is quite 
reasonable for winter and summer circulations. In general, the lower tropo- 
spheric summer circulation features and the upper troposherlc winter circula- 
tion features compare very well with the respective observed circulations. 

This Is not the case for the upper tropospheric summer circulation and the 
lower tropospheric winter circulation features. One of the serious deficien- 
cies of the model both In summer and winter simulations Is the unrealistically 
cold polar upper troposphere and the excessive zonal velocity at the uppermost 
level . 

The analysis of spatial and temporal fluctuations Indicates that, while 
the GCM behaves realistically In many respects, some discrepancies between the 
model and tne atmosphere are evident. With regard to the stationary wave 
field, the zonal variance In the model planetary waves shows a latitude-height 
structure very similar to that observed. However, the GCM's simulation of the 
amplitude and phase of individual planetary waves Is not as realistic. Maps 
of local variability show tliat the model's band-pass variances generally are 
In accord with observations, with the exception of the upper level transient 
heat and momentum fluxes. The low-pass variances are not successfully simu- 
lated however. In the wavenumber-frequency domain, the model Is very skillful 
In handling the medium frequency synoptic scale waves, but has much less 
success In the simulation of the low frequency planetary waves. With regard 
to the latitude-height structure of the total transient fluxes, the heat flux 
Is much more realistically simulated than Is the momentum flux. 

The model's simulation of dlabatlc and frictional processes Is generally 
satisfactory. The simulated evaporation, surface sensible heat flux, and 
surface stress are In good agreement with observations, although the sensible 
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hent flux In somewhat too high. The simulated precipitation distribution is 
realistic and shows the expected predominance of cumulus precipitation in the 
tropics. The most serious deficiency is the excessive frequency of super- 
saturation clouds in the lowest model layer. We believe that the problem is 
due to certain weaknesses in the boundary layer parameterization; appropriate 
revisions are in preparation. 

The simulated Hadley circulations are weaker than observed, but their 
location and depths are satisfactory. The model produces a vigorous subtropi- 
cal Walker circulation in the winter hemisphere. The rising branch of the 
Walker cell coincides with an oceanic precipitation maximum, while the sinking 
branch is over land. 
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